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m ■ ABSTRACT 

Q ■ The kinematics of neutral gas and warm ionized gas have been mapped in one- 

^ ■ dimension across ultraluminous starburst galaxies using the Na I AA5890, 5896 

■^ ■ absorption-line and Ha emission-line profiles, respectively, in Keck II ESI spectra. 

p i' Blue-shifted, interstellar Na I absorption is found along more of the slit than 

Q ■ anticipated, exceeding scales of 15 kpc in several systems. If the outflow diverges 

^— > . from the nuclear starburst, the absorbing gas must reach similarly large heights 

d ■ along the sightline in order to cover so much of the galaxy in projection. If the 

>• . outflow is launched from a much larger region of the disk, however, then the 

S^ . scale height of the absorbing material can be much lower. The large velocity 

gradient measured across a few outflows is, in fact, inconsistent with a flow 
diverging from the central starburst, as angular momentum conservation reduces 
the rotational velocity of an outflow as it expands. In a few young mergers, 
where the galactic nuclei are still separated, the spatial gradient of the outflow 
velocity is shown to be similar to that of the projected orbital motion across the 
system. Widespread star formation triggered by the merger, or the galaxy-galaxy 
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merger itself, may shock heat interstellar gas throughout the disk to generate a 
hot wind, in which cool gas presenting these properties is entrained. The double 
nuclei systems present the highest outflow masses, mainly due to the large area 
over which the cool gas can be detected in absorption against the continuum. 
In a typical ULIG, the mass carried by the cool phase of the outflow is at least 
~ 10^ Mq, which is a few percent or more of the total dynamical mass. Assuming 
the starburst activity has proceeded at the observed rate for the past 10 Myr, 
the kinetic energy of the cool outflows is a few percent of the supernova energy, 
consistent with starbursts powering the outflows. The cool wind is expected 
to be accelerated by momentum deposition, possibly from radiation pressure as 
well as supernovae. The turn-around radii of the cool outflows are at least 30 
to 90 kpc, which presents a significant Na I absorption cross section. If most 
L > O.IL* galaxies pass through a luminous starburst phase, then relics of cool 
outflows will create a significant redshift-path density. Galaxy formation models 
should include this cool phase of the outflow in addition to a hot wind in feedback 
models. 

Subject headings: galaxies: formation — galaxies: evolution — galaxies: fun- 
damental parameters — ISM: kinematics and dynamics ISM: structure ISM: 
evolution 



1. Introduction 

Resonance absorption lines provide a more sensitive probe of diffuse gas in galactic 
halos than do emission lines. The plethora of lines in the rest-frame ultraviolet spectrum 
has made it relatively straightforward to establish the ubiquity of outflows in high redshift 
galaxies (Franx etal. 1997; Pettini etal. 2002; Shapley etal. 2003). Space-based observations 
show these lines are also blueshifted in the UV-brightest nearby galaxies (Heckman & Lei- 
therer 1997, Kunth etal. 1998, Vazquez etal. 2004). Optical resonance lines also probe halo 
gas kinematics and allow measurements in more nearby galaxies; but the paucity of lines 
limits these studies to low ionization states. In particular, interstellar Na I A, A5890, 5896 
absorption is prevalent along sightlines through the Milky Way halo and was recognized 
in spectra of infrared luminous galaxies (LIG's) (Armus, Heckman, & Miley 1989). More 
recently, spectral observations of this doublet have revealed cool outflows in many starbursts 
(Veilleux et al. 1995, Heckman et al. 2000, Rupke et al. 2002; Schwartz & Martin 2004, Martin 
2005; Rupke etal. 2005b; Rupke etal. 2005d). 

One difficulty in the interpretation has been that the most often observed lines trace the 
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kinematics of low-ionization species rather than the hot, i^lO'^ K, gas that powers these mul- 
tiphase outflows. The abundant Na I lines probe regions where H is neutral - i.e. relatively 
cool, dense, and dusty regions. Hydrodynamic instabilities along the wind/disk interface, 
and at the accelerating shell, generate fragments of interstellar gas that are entrained in the 
hot wind. The ram pressure of the hot wind probably accelerates these clouds (Heckman 
etal. 2000, Fujita etal. 2006), although radiation pressure may also be important at ultra- 
luminous luminosities (Murray, Quartaert, & Thompson 2004; Martin 2005; Rupke etal. 
2005b; Rupke etal. 2005d). Martin (2005, hereafter Paper I) showed that some cool gas 
is accelerated to velocities ~ 500 km s~^ in galaxies with SFR > 10 M0 yr^^, but that 
the terminal velocities are lower in dwarf galaxies. A simple dynamical model originally 
suggested by Heckman etal. (2000), and subsequently developed further by Murray etal. 
(2005), was used in Paper I to demonstrate that the cool gas is generally accelerated to the 
terminal velocity of the hot wind by ram pressure; however, some dwarf starbursts do not 
supply enough momentum to bring the entrained mass up to the speed of the hot wind. 
Rupke et al. (2005b) confirm this empirical result but suggest the velocities saturate due 
to a lack of gas remaining in the path of the outflow or a decrease in the efficiency of the 
thermalization of supernova energy. 

This paper presents longslit spectroscopy of the Na I absorption across ultraluminous 
infrared galaxies (ULlGs), defined as galaxies with far-infrared luminosities similar to that 
of quasars - LboI ~ Lm > 10^^ L0 (Sanders et al. 1988). ULIGs may be the local analogs of 
the high-redshift, dusty starbursts that contribute significantly to the cosmic star formation 
rate at z > 1 and dominate the far-IR/submillimeter background (Puget etal. 1996). The 
presence of close nuclei and an overall disturbed morphology indicate ULIGs are merging 
systems. The collision leads to a central starburst and the large-scale shocking of the gas disk 
(Jog & Solomon 1992). The relevance to this study is two-fold. First, the ambient medium 
that the starburst wind encounters is considerably more complicated than that encountered 
by the M82 or NGC 1569 winds because the tidal interaction propels streams of material 
outward and generates regions of infiow (Mihos 1999; Hibbard etal. 2000; McDowell etal. 
2003; Colina etal. 2004). Second, the gas shock-heated by the merger may drive an outfiow 
similar to that driven by supernova-heated gas in the starburst region (Cox et al. 2004; Cox 
etal. 2005). 

Although spatially-extended, interstellar absorption has been detected in a few galax- 
ies previously (Phillips etal. 1993; Heckman etal. 2000), the measurements described here 
provide one of the first systematic studies of the angular extent and velocity gradients of 
the outfiowing, absorbing gas. This spatial information provides a model-independent esti- 
mate of the total mass in the cool outfiow. Velocity gradients in the cool, absorbing gas are 
compared to the kinematics of the ionized gas traced by Ha emission, which helps constrain 
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the location of the absorbing gas along the sightline. Estimates of the mass-loss efficiency, 
momentum deposition, and energy deposition are derived from these data. These results 
provide insight into the appropriate dynamical model for the cool outflows and empirically 
constrain feedback parameters in galaxy evolution models. Due to the complicated, multi- 
phase nature of winds, a dynamical model is required to understand the fate of the material 
in the outflow. 

This paper is organized as follows. Section 2 describes the extraction and measurement 
of the longslit spectra. Section 3 presents the velocity gradients across the starbursts in 
the Na I and Ha lines. The amount of material in the cool outflow is discussed in Section 
4. Implications for the source of the outflow, the evolutionary products of ULIGs, and the 
detectability of cool winds as intervening absorption-line systems are explored in Section 5. 
Section 6 summarizes the conclusions. 



2. Observations 

Echellete spectra of 41 ultraluminous starbursts were obtained 2000 September 19-20 
and 2001 March 26-29 using ESI on Keck II (Sheinis etal. 2002). Galaxies were selected at 
60 /xm from the IRAS 2 Jy sample (Strauss etal. 1992). They have bolometric luminosities 
greater than 5 x 10^^ Lq, a 60;um bump F^{60nin) > F,^(12/im) x Fj,(25/im), and declination 
6 > —35° . The CO-sample contains 15 ULIGs with CO velocities from Solomon et al. 
(1997) and 3 ULIGs with CO velocities from Dr. Aaron Evans (pvt. comm.). The echellete 
has no automatic compensation for atmospheric dispersion and covers a broad bandpass 
(A3900 — l.lyum ), so the position angle of the 20" ESI slit was generally chosen near the 
average parallactic angle. The exceptions were a few objects where structures in R and K 
band images (Murphy et al. 1996) warranted specific position angles; and these systems were 
simply observed at a low airmass. The spectral resolution, about 70 km s~^ , is not high 
enough to resolve individual clouds and shell fragments (Schwartz & Martin 2004; Rupke 
etal. 2005c; Fujitaetal. 2006). 

Paper I describes the extraction, reduction, and analysis of the one-dimensional nuclear 
spectra in the CO subsample. This paper focuses on the analysis of the two-dimensional 
longslit spectra. The reduction was identical to Paper I with the following additions. The 
longslit spectra were rectified independently for each order using arc lamp exposures. The sky 
was fitted along the spatial dimension and subtracted from the co-added data frames. The 
spectral images were examined around Na I and Ha . A series of one-dimensional spectra 
were extracted to enable measurements of the spatial gradients in the Na I line profile. Care 
was taken to account for variations in plate scale among orders when extracting Ha emission 
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from the lower order spectrum of the same physical region. Each order was flux calibrated 
independently using observations of spectrophotometric standard stars (Massey etal. 1988). 



3. Kinematic Results 

The nuclear spectra presented for the CO-observed sample in Paper I reveal cool gaseous 
outflows in 15 of 18 ULIGs. This incidence of winds, ~ 80% of ULIGs at z ~ 0.1, is higher 
than that in local infrared luminous galaxies (Heckman et al. 2000; Rupke et al. 2005b). Since 
the systemic velocities are known better for these ULIGs than others in the 2 Jy sample, 
this CO-sample is also used here to explore the gas kinematics across the merging galaxies. 
The cool component of the outflow is traced in Na I absorption, and the warm photoionized 
component is seen in Ha emission. The absorption can occur at any location along the 
sightline, but the peak of the Ha emission proflle will be dominated by denser gas, which 
likely resides in some form of a disk. Random motions dominate the stellar dynamics in 
ULIGs, but signiflcant rotation is observed in the stellar component; and the gas presents 
even more rotation (Genzel etal. 2001; Tacconi etal. 2002). This decoupling of the gas and 
stellar dynamics occurs because of violent relaxation. By the time the nuclei are 1 kpc apart, 
the stellar systems have nearly reached their equilibrium values of rotation and dispersion, 
but the accretion of high-angular momentum gas continues (Mihos 1999; Bendo & Barnes 
2000). 



3.1. Spatial Extent of Gaseous Outflows 

The Ha emission (and more recently the X-ray emission) from galactic winds have 
provided direct images of bipolar outflows from starburst galaxies (Heckman, Armus, & 
Miley 1990; Martin 1998; Veilleux, Cecil, & Bland-Hawthorn 2005). In constrast to the 
absorption measurements, the faint halo emission is detected beyond the galaxy, usually 
along the minor axis. Indeed, the wind emission is hard to see in projection against the 
galaxy because the disk is much brighter. The dominance of the disk is illustrated by the 
2D-Ha spectra shown in the upper right panels in Figure 1. 

The area covered by the ULIG outflows is explored in a systematic way here. A previous 
(Na I) longslit measurement across the nearby luminous infrared galaxy, NGC 1808, Phillips 
et al. (1993), detected Na I absorption on the near side of the nucleus and scattered emission 
off the back side. The slit spectra presented here are ideal for examining outflow properties 
across the merging galaxies. The nuclei of the younger mergers yield independent high 
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signal-to-noise sightlines, and the stellar continuum emission is sufficiently bright away from 
the nuclei to measure gradients in the projected properties of any outflows. 



3.1.1. Extended Na I Absorption 

Figure 1 shows the Na I spectra extracted at (typically) five positions, which are marked 
on R band images kindly provided by Lee Armus and Tom Murphy. Two ULlGs with 
strong nuclear Na I absorption, IRAS1056+24 and IRAS1720-00, are shown in print; and the 
entire sample appears in the electronic edition. The width of the Doppler-shifted absorption 
component averages < Af^ >= 320±120 km s~^ FWHM toward the nucleus (Martin 2005), 
so the ESI data reveal a great deal of kinematic structure. Many of the off-nuclear sightlines 
present complex line profiles, and the line centers can be measured to ~ 1 km s~^; the 
continuum S/N is adequate over about half of the 20" slit. The absorption-line proffies for 
IRAS1056+24 are a remarkable example of this kinematic structure. Toward the nucleus. 
Aperture 3 in Figure 1, the absorption profile presents a high velocity tail. This component 
is also detected in Apertures 1 and 2 to the west. Surprisingly, the outflow extends beyond 
the main galaxy to the east. It is detected against a companion galaxy in Aperture 5, 7.1 
kpc (8'.'4) from the nucleus. Although the absorption equivalent width declines markedly 
over this distance, the velocity of the line center is remarkably uniform. If we picture the 
absorbing gas residing close to the disk (i.e. within a scale height), it is very difficult to 
understand why the projected bulk flow speed would be this coherent over 11.4 kpc (13'.'4). 
The local star formation rate drops significantly over this distance. Measured U half-light 
radii for ULIGs are typically around 5.3 kpc (Surace etal. 2000). Note also that this is the 
minimum angular extent of the absorbing gas; it cannot be detected at larger distances due 
to the faintness/absence of continuum emission. In the context of an outflow diverging from 
the central starburst region, the absorbing gas would have to be many kpc above the plane 
of the disk given the area covered by it. See §3.2.3 for an alternative interpretation. 

The size of the starburst region has been estimated from the infrared luminosity. The 
tight correlation of lOO/im luminosity and CO luminosity indicates the dust emission from 
ULIGs is optically thick at A < lOO/im (Downes, Solomon, & Radford 1993). A minimum 
dust radius can therefore be estimated from a spherical blackbody at the measured dust 
temperature (Solomon etal. 1997); and the mean size of the starbursts this paper is i? ~ 
200 pc, reaching a maximum of 373 pc in IRAS1945+09. The CO radii listed in Table 2 
are generally larger, presenting a mean of 320 pc and maximum of 575 pc. The starburst 
regions are clearly much smaller than the region where outflowing gas is detected. 

Table 1 summarizes the spatial extent of the absorption. For the 15 galaxies with cool 
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outflows, the Na I doublet is detected over angles ranging from 2" to 13" that subtend 4 to 
18 kpc at the source. Mergers whose constituents are still visibly separated on the sky provide 
some of the widest sightline separations. Toward the object 7'.'3 south of IRAS0018-08, the 
spectrum from Aperture 1 (in Figure 1) presents a marginal pair of Na I lines. Unambiguous 
detections are readily apparent in the spectra of the close companions to IRAS0015+54 (10.4 
kpc), IRAS1150+13 (11.53 kpc), and IRAS1648+54 (5.41 kpc). The continuum S/N drops 
off before the absorption strength weakens in IRAS0803+52, IRAS1049+44, IRAS1056+24, 
IRAS 1524+ 10, IRAS 1720-00, and IRAS2008-03, so the measured extent of the Na I absorp- 
tion is a lower limit for these galaxies. A few systems with extended continuum emission, 
from tidal features, do not present Na I absorption; these are IRAS1836-I-35 Aperture 1 (5'.'7 
S) and 5 (4'.'6 N), IRAS1929-04 Aperture 4 (8'.'1 SE), and IRAS2336+36 Aperture 6 (6'.'6 
NW). 



3.1.2. Extended Ha Emission 

The upper right panels in Figure 1 show the 2D Ha spectra along all the slits. The 
IRAS1056+24 spectrum illustrates some typical features of these spectra. The Ha emission 
comes predominately from one of the two merging galaxies; no Ha emission is detected 
from the companion in Aperture 5. In contrast to the Na I absorption, the Ha emission 
is confined to the center of the main galaxy. The non-detection of extended Ha emission 
from IRAS 1056+24 in consistent with a face-on outflow orientation, which is consistent with 
the morphology in HST NICMOS imagery (Scoville etal. 2000). One expects near face-on 
orientations to favor absorption-line detections of winds, and that edge-on disks with winds 
would be the most likely to show extended Ha and X-ray emission (Lehnert & Heckman 
1996; Strickland etal. 2004a; Strickland etal. 2004b). The gas velocity in the blue wing of 
the Ha line profile is, however, often comparable to, and can exceed, that of the most blue- 
shifted absorbing gas (Rupke etal. 2005b). Ionized gas is likely present in the bulk flows 
traced by Na I absorption, but the Na I and Ha line profiles differ due to the quadratic 
dependence of the latter on the gas density and the requirement that the former be in front 
of the continuum source. 

The spectroscopic slit serendipitously intersected extended. Ha emission in 12 of the 
CO-selected systems. The spatial extent, morphology, and velocities of the extended Ha 
emission are summarized in Table 2. Of particular interest is a Doppler ellipse discovered 
8'.'90 (21.2 kpc) north of IRAS0315+42, which identifies an expanding shell of ionized gas. 
The blue- and red-shifted components of the line profile are separated along the northern 
half of the slit. Using a shell radius of 5'.'88 (14.0 kpc), which is the distance the shell 



protrudes past the bright continuum of the disk, the dynamical age and power requirements 
are large - r ~ 57 Myr and L^jn ~ 2.4 x 10*^^ ergs s""*^ cm'^.-'^ The power requirement is not 
uniquely determined because it depends on the density, n, of the surrounding material that 
is pushed out of the way. A constant star formation rate of 1 M© yr~^ and a Salpeter IMF 
from 1 to 100 M© produces L^ ?5i 7.1 x 10'^^ ergs s~^ in mechanical power from supernova 
(and stellar winds) (Leitherer etal. 1999). Scaling to IRAS0315+42, the maximum power is 
L^ ^ 2.1 X 10"^^ ergs s^^ . Since the average density on these large spatial scales is likely 
< 1 cm~^, this energy injection rate appears sufficient to drive the expanding shell. No shell 
in visible in the continuum image. In the Ha spectrum, it is unclear whether the jet-like 
feature south of the nucleus is the other lobe of a bipolar outflow because the line profile 
has a single peak. These jet-like Ha emission features are also seen northeast of the higher 
velocity component in IRAS1049+44, east of the higher-velocity component in IRAS1609-01, 
and toward the western side of the eastern component in IRAS1150+13. 

The bipolar jets extending from IRAS2336+36 terminate at bright emission blobs sepa- 
rated by 15.4 kpc in projection. Note that velocity shear between the halo blobs is opposite 
to the direction of shear across the disks. The diffuse Ha extending from each end of the 
galaxy pair in IRAS0803+52 and IRAS 1836-1-35 appears to be tidal in origin. They are 
mirror images of each other but originate from respective nuclei rather than a single star- 
burst. The smooth, diffuse halos surrounding IRAS0018-08 and IRAS0026-I-42 could have a 
starburst-outfiow origin or a tidal origin. The 2D spectral images dramatically show that the 
extended emission does not share the disk kinematics in any of these cases. The excitation 
of this extended Ha emission will be discussed in a forthcoming paper. 



3.2. Velocity Gradients 

The Ha emission is weighted toward the higher density gas near the disk because the 
emission measure scales as density squared - EM = J nids. The Ha emission from most 
infrared luminous galaxies and a sizeable fraction of ULIGs present rotation (Rupke et al. 
2005b). Along the slit, the gradient in the Ha velocities at maximum intensity establishes 
a baseline for the disk rotation to which the kinematics of the cool, neutral gas can be com- 
pared. For the CO-subsample, Ha spectra were extracted from the same physical apertures 
as the Na I spectra. The aperture locations are shown in Figure 1 on R band images (Murphy 



^These values were estimated by assuming energy is continuously injected in a homogeneous medium 
(Weaver et al. 1977) . 
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etal. 1996)^ 



3.2.1. Ha Kinematics 

The Ha spectra show that the relative velocity of the merging galaxies helps distinguish 
the two nuclei better than imaging alone. Murphy classified 7 of the 18 galaxies in the CO- 
sample as double nuclei based on imaging with resolution ~ O'.'S (or 1.48 kpc at z=0.1). 
The spatial resolution of the Keck ESI data are similar, but the rotation of the individual 
disks make it possible to separate them in two additional cases. The estimated separations 
are IRAS1150+13 (S'.'O, 11.4 kpc) and IRAS1720-00 (3'.'0, 2.5 kpc). Genzel etal. (2001) 
also classified IRAS 1720-00 as a single nuclei system. Their measurement of the stellar 
kinematics, Vrot/o' = 0.48, is consistent with the mean of their double nuclei control sample, 
however; so the reclassification suggested here is dynamically possible. In general, the nuclear 
separation is a useful indicator of the dynamical age of a merger (Murphy et al. 2001). That 
nine of the 18 ULIGs discussed here present double nuclei indicates half the sample have not 
advanced to the final stage of the merger yet. 

Voigt profiles were fitted to the Ha lines in Figure 1, and Doppler shifts were derived 
from the central wavelength. The Ha position-velocity (P-V) diagrams are shown (by solid 
circles) in Figure 2. The spatial origin is defined as the nucleus of the ultraluminous galaxy. 
Tilted lines reveal projected rotation. Again using IRAS1056+24 for illustration, the Ha 
P-V diagram reveals the ionized gas in the galaxy is receding (approaching) on the western 
(eastern) side of the galaxy. The velocity difference is 120 km s^^ . The low projected 
rotation speed of 60 km s~^ is consistent with a disk seen nearly face-on.^ The nucleus lies 
close to the region where the P-V curve crosses zero velocity. Indeed, the rotation curves 
provide another check on the assigned systemic velocities. The uncertainty in measuring 
the strong Ha line profiles is neglible, but the positions in Figure 1 are a bit fuzzy for 
the following reasons. First, within each aperture, the position of the brightest region is 
not known. A shift perpendicular to the slit could shift the velocity scale by as much as 
38 km s~^ . More importantly, the brightest emission knot could be anywhere in the aperture 
along the slit, so the spatial error bars are half the separation between points. And, finally, 
the heavy extinction toward these nuclei means that the nucleus identified in the R band 
images may be slightly offset from the true nucleus at longer wavelengths. Hence, as long as 



^Exact positions of each aperture can be determined from the slit position angle, see Table 1, and the 
distance from the nucleus, shown in Figure 2. 

■^For example, a galaxy with circular speed of 250 km s^^ inclined 14° will produce the observed gradient. 
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the velocity of the R-band nucleus is close to zero, there is no reason to doubt the systemic 
velocity. 

In the position- velocity diagrams or two-dimensional Ha images, the slope of the emis- 
sion across a single galaxy shows the projected direction of the rotational spin vector. Sim- 
ilarly, the tilt of an imaginary line connecting the centroids of the Ha emission of the two 
galaxies indicates the projected orbital rotation. For example, in Apertures 1, 2, and 3 
for the IRAS1720-00 pair, the Ha velocity shear traces the rotation curve of the northern 
galaxy. Apertures 4 and 5 show the projected rotation of the southern galaxy. Notice that 
the velocity shear between the two galaxies in IRAS1720-00 is of the same order as that 
within the northern galaxy itself. Hence, the Ha PV-diagrams reflect orbital kinematics as 
well as rotational motion. 



3.2.2. Na I Kinematics 

Spectra around the Na I doublet were extracted from the same locations as the Ha 
spectra. These apertures, which are marked in Figure 1, were chosen to optimize sensitivity 
to interstellar absorption. As in Paper I (Martin 2005), the simplest fitting procedure holds 
one pair of lines at the systemic velocity and finds the best-fit velocity for a second pair of 
lines. The absorption troughs in ULIG spectra are more complex than this simple description, 
but this approach has the advantage of characterizing the bulk velocity of the dynamic 
component without introducing a model-dependent components. The fitted velocity depends 
on the doublet ratio, which was initially fitted. As described in §4.0, constraining the 
column densities required a more sophisticated approach; and these spectral fits are shown 
in Figure 1. The velocities measured with both methods are plotted in Figure 2 and made 
available in Table 3. These approaches generally yield consistent velocities, particularly when 
the absorption troughs present significant velocity structure. It was important, however, to 
account fot the systemic absorption component, which was sometimes much stronger than 
that from stellar absorption alone (see Paper I). The Na I fitting was carried out with the 
IRAF distribution of specfit (Kriss etal. 1994). Limitations of the specfit parameter range 
are discussed by Rupke etal. (2005a). In the small number of cases where difficulties were 
encountered in practice, fitting parameters where frozen at a plausible value and fiagged. 
The spatial resolution of the data presented here, typically ~ 1 — 3 kpc, provides several 
velocity measurements across each ULIG, which are used here to explore velocity gradients 
across the outfiowing, absorbing gas. 

The Na I absorbing gas is not found in simple rotation like much of the emission-line 
gas. Across IRAS 1056+24, all the Na I velocities are blue-shifted relative to the systemic 
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velocity. The bulk of the absorption trough is rarely redshifted. One unique exception 
is IRAS0803+52, where all interstellar absorption is redshifted; and the cool gas is either 
infalling or ejected from a foreground component of the system. Section 3.2.3 and Ap- 
pendix A explain how outflows could produce redshifts at a few positions. In systems like 
IRAS1150-I-13, redshifts can also be produced if an extended rotating disk is seen in projec- 
tion against a galaxy behind it. The position - velocity plots in Figure 2 clearly demonstrate, 
however, that the bulk of the absorption trough is generally blueshifted relative to the CO 
(and Ha ) velocities. 

The left panel of Figure 3 compares the velocity gradients measured across all the 
systems in Ha and Na I; these measurements are summarized in Table 2. Many of the 
galaxies, including IRAS1056-I-24, present relatively small Na I velocity gradients, falling 
below the diagonal line in Figure 3. Galactic winds are not expected to exhibit much rotation. 
In most models, outflows diverge from the central region of a galaxy, so angular momentum 
conservation slows down any rotational component present at launch. Appendix A describes 
a simple geometrical model of such an outflow and demonstrates that the projected velocities 
should be nearly constant along a slit aligned with the major axis of the galaxy. Panel a 
of Figure 4 shows several systems with little velocity gradient - IRAS 1524+ 10, IRAS 1929- 
04, and IRAS2336+36. Higher resolution mapping of the gas distribution is required to 
determine whether the slit position angle follows its major axis. 

A significant velocity gradient was discovered across IRAS 1720-00. In the Na I longslit 
spectrum (shown in Figure 1), the gradient is visible in the absorption lines between aper- 
tures 1 and 5. The position-velocity diagram shows that the Na I lines shift by 42 km s~^ 
kpc^^ across the system, which is just as large as the Ha velocity gradient of 47 km s~^ 
kpc^^. Considering the large offset, about 400 km s^^ , between the Na I absorption and 
the Ha emission, the Na I and Ha velocity gradients are remarkably similar. In both Na I 
absorption and Ha emission, the projected orbital motion of the merging galaxies produces 
much of this shear. Tic marks in the position-velocity diagram mark the spatial extent 
of each galaxy along the slit. The velocity gradient between the nuclei of IRAS1720-00 is 
38 km s~^ kpc~^. For all the double nuclei systems that had line detections in both galaxies, 
this orbital shear was measured using the Ha emission line. The right panel of Figure 3 
shows that the total Na I velocity gradient across these is similar to the orbital shear. The 
cool outflow in IRAS0015+54 shows similar kinematics; i.e. the Na I absorption has a large 
blueshift relative to the Ha emission, but both have large velocity gradients that are partly 
due to orbital motion. The absorbing gas cannot be part of the disk in these two systems. A 
rotating gas disk would produce redshifted absorption on one side of the center of rotation. 
The apparent rotation could be explained by invoking separate outflows from each nucleus. 
However, the starburst activity is often much stronger in one nucleus, so it is not obvious 
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that two outflows are produced. 

Some cool outflows present more rotation than expected. Across IRAS1720-00, aper- 
tures 1 and 2 cover the ULIG, and apertures 3 and 4 cross the companion. The P-V diagram 
shows the rotation curve of the ULIG is 65 km s^^ kpc^^, while the companion appears 
to rotate in the opposite direction at 20 km s^^ kpc^^. A similar paradox has been de- 
scribed for halo gas in the spiral galaxies NGC 891 and NGC 5775. The gas velocity falls 
off more slowly with scale height in NGC 891 and (to a lesser extent) in NGC 5775 than 
ballistic models predict (Rand 2000, Tullmann etal. 2000). An outward pressure gradient, 
or drag between the disk and halo, seems to be required (Collins etal. 2002). The situation 
is not completely analogous across IRAS1720-00 and IRAS0015+54 because much of the 
velocity gradient arises from orbital motion. Rotating halo gas is also a common property 
of intermediate redshift {z ~ 0.5) spiral galaxies (Steidel etal. 2002). 

Six systems in the left panel of Figure 3 present a Na I velocity gradient that is larger 
than the Ha gradient. In only two systems, IRAS0018-08 and IRAS1049+44, are the 
velocities not constrained well enough in for the large Na I gradient to be significant. For pre- 
mergers, the velocity gradient due to the projected rotation of the galactic gas was measured 
from Ha velocity gradient across the ULIG member of the galaxy pair. The middle panel of 
Figure 3 compares the total Na I velocity gradient to this rotational gradient. The similar loci 
of the data in the middle and left panels indicates that much of the Na I velocity gradient 
reflects rotation of the galactic gas. When the slit is aligned along the minor axis of the 
outflow model from Appendix A, high inclination observations produce a velocity gradient 
across the center of the galaxy. Four systems are shown in panel b that could plausibly reflect 
this type of break in their P-V diagrams. The absorption velocities become redshifted when 
the outflow cone is pointed away from the observer on the far side of the minor axis. For 
the case of IRASl 150+13, an alternative explanation is that the receeding western side of 
the gas disk associated with the eastern galaxy is being detected in absorption against the 
western galaxy. This geometry naturally explains the redshift, but the Na I column densities 
do not decrease as rapidly as might be expected for a disk seen in projection. Another origin 
is needed for the rotation in IRAS0315-I-42 and IRAS1836-f35 because the Na I velocity 
gradient is in the opposite sense of that of the star-forming disk.^ 

Spatially-extended starburst activity throughout both gas disks would explain the oberved 
rotation of the outflow. If the the wind is generated over a large fraction of the disk, then 



''The gradient is unlikely minor axis shear. The IRAS0315+42 system is thought to be viewed edge-on 
in light of the detection of the Ha shell, but it would need to be viewed at low, i < 90° — 6*^, inclination to 
keep the velocities on both sides of the nucleus blueshifted as observed. 
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the absorbing gas need not be at extremely large heights to completely cover the continuum 
source. Higher resolution studies of the star formation activity are needed to confirm the 
viability of this explanation, however. The size of the molecular gas disk and the far-infrared 
emission are much smaller than the area covered by the outflows described in this paper. 



4. Mass of Cool Gas in Outflow 

Measurements of mass outflow rates are important for understanding the influence of 
winds on galactic evolution and the intergalactic medium. The mass of hot and warm gas in 
winds has been shown to increase with the SFR (Martin 1999). The total mass in the cool 
component is difficult to estimate. The strength of the Na I absorption along the slit, and 
its relation to the amount of cool gas along the sightline, are discussed here. 



4.1. Na I Column Densities 

Extraction of the Na I column density is complicated by the blending of the AA5890, 5896 
lines and their saturation. The inferred columns grow non-linearly with the line strength in 
this regime, so seemingly minor degeneracies in profile fitting can lead to enormous allowed 
ranges for the estimated columns. A comparison of standard curve of growth techniques 
(Spitzer 1968) and more specialized methods was therefore required. For a related but 
distinct approach to fitting absorption troughs in ULIG spectra, see Rupke et al. 2005a. 

The large width of the Na I absorption trough in ULIGs indicates the presence of mul- 
tiple absorbing regions along the sightline; see, for example, the simulations of Fujita etal. 
(2006). The underlying velocity distribution is likely more complex than the Gauassian 
profile assumed in the standard curve-of-growth method. Simple deconvolution techniques 
were therefore intially applied to the line profiles using a kernel composed of two narrow 
lines of fixed intensity ratio, which are separated by 300 km s^^ . The resulting templates 
had Gaussian-like features including, for example: (1) A single Gaussian centered at the 
systemic velocity, (2) A smooth distribution with a maximum at the systemic velocity but a 
strong skew toward blueshifted velocities, and (3) two Gaussian profiles - one at the systemic 
velocity and one Doppler shifted. The ambiguity in the kernel can, in principle, be elimi- 
nated by working on the apparent optical depth profile instead of the intensity profile, as 
T589o(A) = T5896(A + s)(/A)589o/(/A)5896- Substituting for the separation of the doublet s and 
the oscillator strengths, this expression becomes T589o(A) = 1. 98x5396 (A(A) + 6.97(1 + z)). In 
practice, the optical depth is only a simple function of the intensity profile when overlapping 
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absorbers completely cover the continuum source. The constraint imposed by the doublet 
can also introduce oscillations in the solution (Arav et al. 1999a), and regularization methods 
have not been developed for the case of partial covering. 

Since the templates suggested by this exercise are well-behaved, pairs of Gaussian pro- 
files were fit to the spectra following Heckman etal. (2000), Rupke et al. (2002), and Martin 
(2005). When the underlying template can be described by many Gaussian components, 
and the invidual components are well-behaved - i.e. Tq< 5 — 10 and the Doppler parameters 
are not varying wildly, multiple components can be analyzed collectively using standard, 
single-component curve-of-growth methods (Jenkins 1986). In Paper I, one component was 
required to be at the systemic velocity; and outflow/inflow was associated with one (or more) 
additional components as required to described the overall profile shape. This approach con- 
strained the doublet ratio well provided the line profiles had multiple absorption troughs; 
but the doublet ratio had to be guessed for systems presenting smooth absorption troughs. 
The doublet ratios listed in Table 3 of Paper I imply very high optical depths, tq > 10, for 
6 systems - IRAS0015+54, IRAS0018-08, IRAS 1049+44, IRAS1056+24, IRAS1720-00, and 
IRAS1836+35. (Here, tq is the estimated optical depth at the center of the A5896 line.) 

The curve-of-growth method clearly breaks down when the lines are strongly saturated. 
The accuracy of alternative methods, like integrating the apparent optical depth, are also 
limited by the degree of saturation. The spectral resolution of 70 km s~^ may not resolve 
individual components, and the only way to ensure accurate measurements is to obtain much 
higher resolution spectra in the hope of resolving out the saturated regions. For now, one 
self-consistent approach is to measure a lower limit for the optical depth and then fix the 
doublet ratio to give this same optical depth. 

In practice, extraction of the apparent optical depth is complicated by the non-unity 
covering fraction of the absorbing gas. High-resolution, ultraviolet spectra of both high- 
redshift and nearby starbursts often present completely black, low-ionization lines (Pettini 
etal. 2002; Vazquez etal. 2004). Because the doublet ratio indicated the Na I lines are 
optically thick, the non-zero residual intensity can be attributed to partial covering of the 
continuum source. A single aperture averages sightlines over areas of a square kiloparsec, 
which are likely much larger than the dense regions where Na I is abundant. Along sightlines 
through our Galaxy, the Na I doublet reveals small scale structure, around 3200 AU (e.g. 
Andrews, Meyer, & Lauroesch 2001). The difference in covering factor between the optical 
and UV lines may simply reflect the patchiness of the UV light and its association with 
outflow regions. When the doubet lines are not blended, the covering factor can be derived 
directly from their residual intensities (Barlow & Sargent 1995; Hamann etal. 1997). To 
extract the optical depth from the intensity profile, a velocity-independent covering factor 
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was simply estimated from the residual line intensity as C/ ~ 1 — min{I{X)). The aver- 
age value for the nuclear apertures are Cf ~ 30% (Martin 2005). If the overlapping line 
components are ignored, and the /5890 measured from the bottom of the observed profile, 
the mean increases to Cf ~ 45%; these values are used here for comparison to Heckman 
etal. (2000). Galaxies with larger covering factors nearly always show stronger Na I lines 
(Heckman etal. 2000; Martin 2005). Inspection of Figure 1 shows that the nuclear covering 
factor varies among galaxies, and that the covering factor typically decreases with distance 
from the nucleus along any particular slit. 

An apparent optical depth profile, r(A), was extracted from each spectrum. When 
multiple, overlapping absorbers are present, the intensity of the absorption trough, /(A)//o = 
Cfe~'^ + 1 — C/, is determined by the amount of light not covered by absorbers, (1 — Cf)lQ, 
and the optical depth in each transition at any wavelength, r(A) = r589o(A) + T5896(A). 
The integrated optical depths obtained with this procedure lie in a fairly narrow range, 
2 < T5895 < 10. Because upper limits were adopted for the absorber covering fraction, 
the estimated columns may be systematically low. This decomposition is the simplest one 
consistent with the data. See Rupke et al. (2005a) for a description of other possible absorber 
geometries. 

The optical depth of the Na I A5896 transition was used to estimate the column density 
of the ion. The apparent column density, 

rrieC^ 
or 



Nap = ^VT2 / ^A^A, (1) 



Nap{NaI) = 1.02 X 10^^ cm^^ / r5896(A)rfA (2) 



for the A5896 transition, is a lower limit on the true column. The results are given in 
column 6 of Table 3. They generally agree, to within a factor of a few, with the column 
densities obtained by fitting the intensity profile, as shown in column 8 of Table 3 . Four 
of the ULIGs are included in the low-redshift Rupke et al. (2005a) sample. The estimated 
columns for the nuclear aperture are within a factor of two of their values, which is very 
good agreement considering the different fitting techniques, possible aperture differences, 
and nonlinear response of the columns. As a cautionary note, columns 11 and 12 show the 
minimum and maximum column densities consistent with the intensity profiles. These ex- 
trema were obtained by fitting the intensity profile and allowing the doublet ratio to vary, 
and they illustrate the potentially large uncertainties that arise from the non-linear relation- 
ship between the line profile and the optical depth. In addition, it should be recognized that 
the absorption troughs represent averages across the continuum source. Since the non-linear 
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weighting strongly favors low optical depth paths, even techniques that fit the optical depth 
will tend to underestimate the column density. 



4.2. Estimated Mass in the Cool Wind 

Estimating the H gas column N{H) from the Na I column requires corrections for 
the abundance of Na relative to H, the depletion of Na onto dust grains, and the amount 
of ionized Na. Sub-solar abundances measured along sightlines through the disk and halo 
of the Milky Way have been interpreted as depletion onto dust (Wakker & Mathis 2000; 
Savage & Sembach 1996); and a solar Na/H abundance ratio, 2.04 x 10~^ (Morton 1991 and 
references therein), is adopted for these massive ULlGs. A conservative conversion to total 
HI column, 

20 __2 f NjNal ) 
W^ cm 

UN{Na)/N{NaI)) 



N{HI) = 4.90 X 10^° cm-2 : \ x 



10 



(3) 



is applied here, where fd is the depletion factor d^a = dx = {X/H)q/{X/ H). From the 
range of measured N(Na I) in Table 3, the corresponding H columns are estimated to range 
from N[H) ^ 8.8 x 10^^ cm~^ in IRAS1056+24 (aperture 4 component 1) to as much 
as N{H) ^ 1.2 X 10^^ cm~^ in IRAS0026+42 (aperture 1). In many regions where H is 
neutral, Nal is probably not the dominant ionization stage. (The ionization potential of 
sodium, 5.1 eV, is much less than that of HI.) In the Milky Way, the product of the Na 
depletion onto dust grains and the Na I ionization fraction presents little variation with 
N(H). Provided the H column is less than 10^^ cm"^ , their product is ~ 300, within a 
factor of a few (Wakker & Mathis 2000). Since the correction may decrease by an order of 
magnitude at higher columns due to self-shielding of Na from ionizing radiation (Ferlet et al. 
1985), a correction factor dN idion = 10 is adopted here for ULIGs to avoid overestimating 
the mass column. Using the mean values for intermediate velocity clouds in the Milky Way 
would raise all masses quoted here by a factor of 30. Note that Schwartz & Martin adopted 
dNaidion = 100 in their study of dwarf starburst galaxies because they have significantly 
lower Na I columns. 

The run of column density along each slit constrains the total mass of the cool outflow 
better than the nuclear sightline alone. To estimate the outflow mass, the column density 
at a particular position was multiplied by the area of an annulus, centered on the nucleus, 
and the estimated gas covering fraction. The radius was set equal to the aperture's distance 
from the nucleus, and the thickness of the annulus determined by the aperture spacing. The 
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masses from concentric annuli were added to arrive at the total mass of cool gas in the 
outflow. A significant fraction of the total mass is detected in the off nuclear sightlines. In 
several galaxies, the absorption is strong out to the limit of the continuum detection, so the 
areas are lower limits; and masses are given as lower limits in Table 1. 

Figure 5 shows the mass of the outflow derived this way. Nearly all the error ranges are 
consistent with M,. ~ 6 x 10^ Mq in the outflow. Although no simple relation with the star 
formation rate is observed, the baseline in SFR is not large; and the lack of correlation is 
not surprising. Similar conclusions about the cool wind mass in ULIGs have been reached 
by fitting free-streaming outflow models (constant velocity and mass flux) to nuclear spectra 
(Rupke et al. 2005b). The mass estimates made here contain conservative assumptions about 
line saturation and the Na ionization/depletion correction. It is therefore plausible that the 
masses are actually systematically higher than shown. In addition, if the mass is calculated 
using the the extreme limits on the column density, or fdfion ~ 100, then values of ~ 10^ M© 
are allowed in many ULIGs. 

Figure 6 shows that the cool outflows carry a signiflcant amount of mass. The molecular 
gas mass ranges from 4 x 10^ Mq to 15 x 10^ Mq for these galaxies (Solomon et al. 1997), and 
the interstellar gas is predominantly molecular in the ULIG phase. The mass inferred for the 
cool outflows are typically a few percent of the dynamical mass in the starburst region, where 
Rco ~ fc-w X 100 pc. The Na I column densities measured toward dwarf starburst galaxies 
are lower than those in ULIGs, and the outflow speeds are also lower (Schwartz & Martin 
2004). Using the atomic gas mass referenced by Schwartz & Martin, and adding molecular 
gas masses in the central regions (Hunter etal. 1996, Scoville etal. 1991, and Weif5 etal. 
2001) yields the points for dwarf starbursts shown by the open symbols in Figure 6. The 
dwarfs present lower columns of blueshifted Na I. Even when larger ionization corrections 
are assumed, e.g. fdfion ~ 100 in Schwartz & Martin (2004), the mass of outflowing gas is 
considerably smaller in dwarf starbursts than it is in ULIGs. The fraction of the interstellar 
gas mass in the cool wind is similar to that measured in the ULIGs. The interstellar gas 
masses for the dwarfs may be overestimated in that some of the HI gas is not associated 
with the starburst region, and this correction would shift the outflows in dwarf starbursts 
to higher gas fraction. The galaxy that appears to lie at 100% efficiency is NGC 4449; and 
the wind mass is likely an overestimate as this component was the most saturated in the 
Schwartz & Martin (2004) sample. 
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4-2.1. Dust Associated with the Cool Gas 

When optically thick gas only partially covers the continuum source, the Na I line 
strength is sensitive to the velocity spread of the absorbers and their covering factor (in 
addition to the optical depth). Previous work has shown that the interstellar Na I strength in 
starbursts is more correlated with covering factor than with line width (Heckman et al. 2000; 
Martin 2005). Since the reddening also tends to be larger toward luminous infrared galaxies 
with strong Na I absorption lines (Armus, Heckman, and Miley 1989, Veilleux etal. 1995), 
Heckman et al. (2000) suggested that more reddened sightlines have larger covering factors 
of cool gas and dust. For the ULIGs discussed in this paper, the reddening was estimated 
using the Balmer decrement; and Table 2 lists the results. The measured H/9 /H7 flux ratio 
consistently indicated higher reddening than did the Ha /H/3 ratio; but the discrepancy is 
attriubted to underestimating the strength of the broad, stellar absorption at H7 , which 
was fit simulataneously with the emission component. Figure 7 confirms that the reddening 
tends to be larger in ULIGs with strong Na I absorption. Of the 5 highest equivalent width 
ULIGs, the covering factor is high in four systems - IRAS1056+24, IRAS1524+10, 11720-00, 
and IRAS1836+35. The absorption in the other galaxy, IRAS2008-03, presents an average 
covering factor but the largest velocity width in the entire sample. Figure 8 shows that all 
the galaxies with high covering factor are indeed highly reddened. However, some highly 
reddened starbursts have average Na I covering factors. The ULIGs and LIGs present the 
same range of covering factors, but the ULIGs have higher average reddening. 

It is interesting to examine how much of the reddening might be associated with dust 
in the wind. Along sightlines through the solar neighborhood, the amount of reddening is 
correlated with the gas column (Spitzer, Figurel.l; Savage etal. 1977). Extrapolating this 
Galactic gas-to-dust ratio to the outflowing gas columns in ULIGs yields considerably less 
reddening than is measured. Even for the sightlines with the largest columns, N{Nar) ~ 
10^^ cm~^ , the correction factor fdfion would have to be much larger than 10 for all the dust 
to be associated with the the outflowing gas column. It is not surprising that the reddening 
is larger than that associated with the outflow because there is a signifiant gas column at the 
systemic velocity; and this is consistent with the variation in Na I covering factors among 
highly reddened ULIGs and perhaps the offset between LIRGs and ULIGs. 



4-2.2. Relation of Wind to Merger Process 

The estimated outflow masses are more uncertain than the column densities because 
the area covered by the absorbing gas is not fully constrained. The Na I doublet can be 
detected at larger radii in dynamically younger mergers because these mergers have wider 
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separations (larger continuum baseline). Four of the five galaxies with mass estimates above 
10® Mq are double nuclei systems. Six of the eight double nuclei systems, IRAS0015+54, 
IRAS0018-08, IRAS1150+13, IRAS1056+24, IRAS1720-00, and IRAS2008-03, have large 
estimates for the outflow mass, upper limits of ~ 10^ Mq. (For the other two double nuclei 
systems, IRAS0352+00 and IRAS 1648+54, no Na I absorption was measured toward the 
nucleus.) The extent of the absorbing gas is at least 5 to 18 kpc in the double nuclei 
systems, whereas absorption is detected in the later stage mergers to a radius of only 4 
to 7 kpc. The single nuclei systems, particularly those without obvious tidal features - 
IRAS0315+42, IRAS0803+52, and IRAS1929-04 in particular, do appear to have lower mass 
outflows. However, the difference in estimated mass arises entirely from the apparently 
smaller area covered by the outflow, as the column densities toward the nuclei do not differ 
appreciably (in moderate spectral resolution at least). Rupke etal. (2005b) also find no 
statistically difference among outfiows in pre-mergers, mergers, and old mergers. 

The separation of galactic nuclei and tidal features are good indicators of the relative 
dynamical ages of ULIGs (Murphy etal. 2001). Outflows are expected to evolve with the 
progression of the merger. Galaxy-galaxy mergers likely trigger multiple phases of star for- 
mation, starting with a nuclear starburst but later progressing across the disk (J. Barnes, pvt 
comm). Nuclear activity is also expected to be exposed as these mergers advance (Hopkins 
etal. 2005). Measuring the evolution in the outflow column density is complicated, how- 
ever, by the sensitivity of the sodium ionization corrections to the hardness of the radiation 
fleld (Murray etal. 2006). Although the only classified Seyfert nucleus in the CO-subsample 
is IRAS1150+13, several of these ULIGs have LINER-like optical line ratios. Overall, the 
column density measurements are not precise enough to explore such variation. However, 
ULIGs with clear signs of nuclear activity do show systematically higher outfiow velocities 
(Martin 2005; Rupke 2005d). Also, ULIGs with starburst /AGN-composite spectra show 
some evidence for a contribution to the momentum from radiation pressure (Rupke et al. 
2005d, Fig. 12b discussion). 



5. Discussion 

Galactic winds are ubiquitous in starburst galaxies, and the mechanical power supplied 
by temporally (and spatially) correlated supernovae is sufficient to power them (Heckman 
etal. 1990). The empirical correlation between the mass loss rate of warm/ hot gas and the 
star formation rate supports this scenario (Martin 1999). The detection of neutral gas in 
starburst outfiows came as a surprise (Heckman etal. 2000; Rupke etal. 2002; Schwartz 
& Martin 2004; Martin 2005; Rupke etal. 2005b; Rupke etal. 2005d), because it is not 
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obvious that such dense clouds survive in a wind. This paper takes a closer look at the cool 
component of outflows in ULIGs. In § 5.2, it is argued that the mass loss rates are signiflcant 
but may not be as large a fraction of the SFR as was estimated for the warm and hot phases 
(Martin 1999), a point made previously by Rupke et al. (2002, 2005b) using a complementary 
approach.^ This conclusion is preliminary for two reasons however. First, the wind efficiency 
in the hot and warm phases is not empirically established for ultraluminous galaxies. Second, 
the mass of the cool outflows could be greatly underestimated if some unresolved regions 
of the absorption troughs are highly saturated. While the data are largely consistent with 
supernova feedback accelerating the cool outflows via the ram pressure of a hot wind, several 
alternative acceleration mechanisms are not yet ruled out. In particular, the viability of 
merger-induced winds is explored in § 5.1 in light of the rotation observed across several 
outflows. 



5.1. Supernova-Driven Winds or Merger-Induced Winds? 

Recent theoretical work suggests outflows might be driven directly by galaxy mergers 
(Cox etal. 2004). The collision of gas-rich galaxies may generate strong shocks, which heat 
the interstellar gas at the expense of orbital energy. The high pressure would drive an outflow 
in much the same manner as the hot interior of a superbubble does. These merger-driven 
winds are not easy to distinguish from supernova-driven winds because the shock velocities 
are few x 100 km s"^ in both cases, generating hot X-ray emitting gas at 10^ to 10^ K. 

The primary distinguishing factors are (1) the spatial scale of the energy-injection region 
and (2) the relation between hot gas mass and the star formation rate. The size of the 
energy injection region will be larger for merger-driven winds than for winds driven by a 
nuclear starburst. The latter develop a bipolar shape because the higher density of the 
disk stalls the shock in the disk plane (de Young & Heckman 1994). If the energy source, 
however, is comparable to the size of the disk, the outflow is more spherical (Cox et al. 2005). 
The morphology of X-ray surface brightness contours may eventually help distinguish these 
geometries. The high fraction of ULIGs presenting outflows in Na I, 16 of 18 objects in the 
CO-subsample and a similar fraction in Rupke etal. (2005b), requires a wide opening angle 
for the outflows; but the large angular extent of the blueshifted absorption regions indicates 
a cylindrical or spherical, rather than conical, outflow geometry. The rotation detected in 



^Notc that Rupke etal. 2005b quote SFR's for the mass range 0.1 to 100 Mq but model energy and 
momentum injection using Starburst 99 models, which normalize the SFR to the mass range from 1 to 
100 Mq. These SFR scales differ by a factor of 2.55 complicating comparisons to the work presented here. 
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a few cool outflows, IRAS1720-00 and IRAS0015+54, would also be consistent with a wind 
geometry similar to a very thick disk. These geometrical constaints could well reflect an 
extended, global starburst rather than a merger-driven wind. If merger-driven winds are 
important at all, energy arguments suggest it is during mergers of massive, gas-rich galaxies 
like ULIGs. Cox et al. (2004) found that the temperature and amount of shocked gas varied 
strongly with the type of orbit and impact parameter, but the total gas consumption was 
only weakly correlated with the orbit. Hence, the properties of the hot wind will vary greatly 
among ULIGs with similar SFRs if the bulk of the heating is due to merger shocks. 

The measured velocities and masses were combined to describe the kinetic energy of 
the bulk flow as KE ^ 0.5^2 Mivf, where the sum was over the circular annuli. This 
observed kinetic energy is strictly a lower limit since the velocities are only line-of-sight 
components. Table 1 summarizes the results. For IRAS1056-f24 and IRAS1720-00, for 
example, the measured kinetic energy is ~ 5 x 10^^ ergs and ~ 4 x 10^^ ergs, respectively. For 
comparison, a Starburst 99 (Leitherer etal. 1999) population synthesis model (Salpeter IMF 
and solar metallicity) normalized to the star formation rate of IRAS1056-I-24 and IRAS1720- 
00 indicates that over a period of 10 Myr, supernovae and stellarwinds supply a mechanical 
energy of 3.0 x 10^'^ ergs and 7.6 x 10^^ ergs. Figure 9 confirms that cool outflows typically 
carry a few percent of the supernova energy. The starburst phase may last longer, perhaps 
100 Myr as favored by some dynamical models (Murphy et al. 1996), increasing the available 
feedback energy and proportionately reducing the inferred efficiency of converting it into 
bulk motion of the cool wind. No evidence for a higher thermalization efficiency in ULIGs 
relative to dwarf starbursts, as suggested by Rupke etal. (2005b), is seen in Figure 9. This 
discrepancy is addressed further in §5.2.1 as the two studies agree on the relevant conlusion 
for this section. Adiabatic, wind-driven bubble models can transfer close to 20% of the 
supernova energy to the kinetic energy of the swept-up shell (Weaver etal. 1977); hence 
starbursts are more than capable of powering the observed outflows. Although the orbital 
energy of the merger is larger than the measured kinetic energy of these cool outflows, the 
thermalization efficiency of merger-driven winds has yet to be established. Considering that 
the supernova-driven wind model successfully explains many properties of the warm and 
hot outflows observed in starburst galaxies, it remain the favored power source for the cool 
outflows, although the momentum supplied may be more important than the energy as 
described below. 
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5.2. Gas Dynamics of Cool Outflows 

How far galactic outlows travel is of great interest. In addition to regulating the galactic 
gas mass, galactic winds are thought to pollute the intergalactic medium with heavy elements. 
In Paper 1, the terminal velocities of the cool outflow were shown to be near the escape speed. 
On the one hand, little gas is detected that is moving fast enough to escape the halo potential. 
On the other hand, once material is moving fast enough to leave the galaxy, it will not be 
further accelerated. Understanding the fate of the outflowing material therefore requires a 
dynamical model. Empirical data for cool outflows should help us identify the best model. 

Murray et al. (2005) show that measurements of mass loss efficiency would, in principle, 
distinguish winds driven by momentum deposition from energy-driven winds. If a fraction ^ 
of the mechanical energy supplied by supernovae and stellar winds is thermalized, then the 
maximum mass loss rate in an energy conserving wind is O.SM^f^ = C,Lw, where Voo is the 
terminal velocity of the wind. Taking C, = 10%, the mass flux in an energy-conserving wind 
relative to the SFR is 

M,„ /300kms~^\^ 

0.985^0.1 , (4) 






oo 



where the mechanical luminosity was taken from a population synthesis model with continu- 
ous star formation for a period of 10 Myr, a Salpeter initial mass function, and upper/lower 
mass limits of 1 Mq and 100 M0(Leitherer etal. 1999). The coefficient would be 2.7 times 
larger at an age of 100 Myr. The cooling time of the gas associated with the Na I absorption 
is likely shorter than the dynamical timescale of the outflow. 

Cool clouds of gas and dust are thought to be accelerated by the ram pressure of 
the hot wind and/or radiation pressure. For a total momentum deposition rate psB, the 
maximum mass flux is Mc, where McV^o = ^Psb- Using the mass flux and energy deposition 
rates from Starburst 99 (Leitherer etal. 1999), the estimated momentum deposition rate is 
PsB = v'2LwM, which is 2.1 xlO^'^ dyne at an age of 10, or about 2.4 times higher at an 
age of 100 Myr. Since some of the supernova and stellar wind momentum vectors will cancel 
out in the starburst region, this estimate is an upper limit on momentum production. For 
a burst duration of 10 Myr, the normalized mass outflow rate of the momentum-conserving 
solution is 



Mr /300 km s 



-1 



Lie . (5) 

M, \ Voo J 

In principle, it should be possible to determine whether the cool wind is momentum-driven 
by measuring the ratio of the mass loss rate to the SFR, i.e. the mass loss efficiency rj, as 
a function of escape velocity (or terminal velocity). These energy-driven and momentum- 
driven wind models are illustrated for reference in Figure 10 for comparison to the data. 
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5.2.1. Outflow Rates in Cool Phase 



The mass-loss rate is Mc{r) = Qfi'mHnH{r)vr'^ for a constant velocity, radial outflow 
into solid angle Q. For a mass flux independent of radius, the column density through the 
wind is 






' ' max 



Assuming a spherical wind, and taking /i = 1.4 to include the helium mass, the mass flux in 
the cool wind is 

M = 141 Me yr- ( JL) (_IJL_) ^ 



400 km s~V V47r 

Column 7 of Table 1 shows the mass loss rates obtained using the column densities from 
col. 6 of Table 3 and the size of the starburst region from col. 3 of Table 2. The mass loss 
rates (in the cool phase) reach 10% of the star formation rate in many ULIGs.^ The mass 
obtained by integrating this wind density over an outflow of opening angle Q from a launch 
radius tq to r^nax is 



M, = 6.89 X 10^ M^ ( -^] (r_rna^c^\ ( ^\ ^ 



1 kpcy \ 1 kpc / v47ry 

4.9 X 1020 cm- 
Taking the starburst radius from col. 3 of Table 2 and the maximum radius from col. 4 of 
Table 1, which is a lower limit on the extent, yields the mass estimates in col. 8 of Table 1. 
The results are similar to the masses obtained by direct integration of the column density. 

These mass loss efficiencies are normalized by the star formation rate (of 1 to 100 Mq 
stars) and plotted with the analogous values for the Schwartz & Martin (2004) galaxies 
in Figure 10. The terminal velocities, i.e. the most blueshifted part of the absorption 
trough, are much lower in the lower mass galaxies than in the ULIGs. Both sets of data 
are consistent with the tracks for momentum-driven winds with e ~ 10 to 100%, but the 



^This statement assumes a lower mass cut-off of 1 Mq for the initial stellar mass function and 
dln{N) / dln{M) = -2.35. 
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lower efficiencies appear to be favored in the ULIGs. Non-unity values for e could reflect 
cancellation of the momentum vectors from multiple supernovae in the starburst region or 
that the hot phase of the wind carries more momentum. The variation of e with mass is 
not well-determined. In addition to uncertainties about the sizes of the absorbing shells in 
dwarfs, the ionization/depletion corrections are highly uncertain. In particular, Schwartz 
& Martin adopted fdfion = 100, but a value of 10 is used here for the ULIGs. These 
assumptions yield, for purposes of illustration, a higher mass-loss efficiency, rj = Mc/M^,, in 
the dwarf starbursts than for the ULIGs in Figure 10. Constant fdfion, in contrast, indicates 
similar momentum deposition efficiency, e ~ 10%, in all the starburst galaxies and reduces 
the variation in mass loss efficiency. 

Rupke et al. (2005b) agree with these results for the efficiencies of momentum, energy, 
and mass deposition in the cool phase of outflows (in the order of magnitude sense or better) 
. However, their paper also claims more efficient energy deposition in larger starbursts. In 
Figure 10, the higher value of fdfion = 100 for the dwarf galaxies yields constant thermal- 
ization efficiencies .^ ~ 1%. Lowering the ionization and depletion corrections in the dwarf 
starbursts (to match the ULIG values) qualitatively reproduces the trend noted by Rupke 
etal. (2005b). The scalings cannot be compared quantitatively because Rupke etal. do not 
consistently use the same lower mass cut-off for the ULIG SFR's, the dwarf galaxy SFR's, 
and the star formation rate of the population synthesis model.. The robust results from 
Figures 6, 9, and 10 and Rupke et al. (2005b) are the efficiencies of mass, energy, and mo- 
mentum deposition in the cool phase of ULIG outflow. More measurements are needed for 
the dwarf starbursts, although it appears at least plausible that the efficiencies do not vary 
much. 



5.2.2. Turn Around Radii 

To explore how far the cool clouds in ULIG outflows might travel, the momentum- 
driven-wind models described by Murray et al. (2005) were used to estimate turnaround 
radii. The momentum deposition falls off as r^^ for both ram pressure driven winds and 
winds driven by radiation pressure. Maximum acceleration takes place near the base of the 
free-flowing wind, and the gravitational potential of a roughly isothermal halo will decelerate 
the flow at large radii. 

Consider the motion in an isothermal halo of circular velocity Vc = 250 km s~^ , or 
velocity dispersion a = 177 km s~^ . For ram pressure driven winds, the clouds begin to 
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decelerate at a radius, 



M, /177kms"^'^ 



R„ = 37 kpc X (9) 

^ ^ 100 M0 yr-i V ^ / 

2.5 X 10^0 cm-2 ^ ^ 
N-„ ■ (1°) 

Integrating the equation of motion from a launch radius of 1 kpc, the turnaround radius is 
found to be so large that the clouds are not likely to fall back. For radiation pressure driven 
winds in the optically-thin limit, the deceleration radius is 

L /177kms"^\^ 



fi» = 3.6 Upc^^^ ^^^ j X (11) 

650 cm2 g-i ■ ^ ^ 

Integrating the equation of motion, one finds the turnaround radius is about 32 kpc. 

It is not currently known which of the above wind models, or an optically thick radiatively- 
driven wind model, is a better description of cool outflows in ULIGs. However, they all have 
large enough turn-around radii to present a significant geometrical cross section for absorp- 
tion. To illustrate the space density of winds, consider a fiducial turn-around radius of 30 kpc 
for the cool outflow, which is a conservative, lower limit based on the dynamical model above. 



5.2.3. Intervening Absorption from Cool Winds 

The association of intervening absorbers with galaxies is well-established for strong, low- 
ionization (e.g. Mg II) systems. The open question is whether these absorbers are accreting 
halo gas or starburst winds. The properties of low-ionization gas in outflows was not easy to 
predict, however, due to the non-equilibrium ionization state in winds and the multi-phase 
fluid. Although the observations presented in this paper are taken along sightlines with a 
special orientation, they do provide some constraints. 

In the local universe, the space density of infrared bright galaxies exceeds that of other 
classes of galaxies for LboI > 10^^ L© (Soifer etal. 1986; Soifer etal. 1987). Adjusting the 
Soifer etal. luminosity function to Hq = 70 km s~^ Mpc~^, the local number density of 
galaxies brighter than 10^^ Lq is 1.0 x 10^^ Mpc^^. For an absorber radius of 30 kpc, the 
redshift path density is small, dN/dz k, 1.4 x 10"'^, and approximately equal to UQacrosscHQ^ 
(within 15% at z=0.1). If the infrared active phase is, however, a nonrecurring event, the 
fraction of normal galaxies passing through this stage approaches unity for infrared-bright 
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periods lasting less than 100 Myr (Soifer etal. 1987). For cool ULIG winds to produce a 
significant number of intervening absorbers, the winds must hang around the galaxy long 
after the starburst ends. 

These relics should be observable as intervening Mgll systems. Steidel & Sargent (1992) 
found the redshift density of strong, Wr > I.OA, Mg II systems was dN/dz = 0.27 over 0.2 < 
z < 2.2. A source density of ng = 0.022 Mpc~^(30 kpc/r)^ is required to create a comparable 
path density. For the space density of normal galaxies to be this high (Blanton etal. 2003), 
galaxies as faint as O.OllL* would need to contribute to the relict wind population. This 
estimate can be improved by applying the measured evolutionary dependence of strong 
Mgll systems (Prochter etal. 2004). Taking into account the increase in path density with 
redshift, the extrapolation to redshift 0.1 reduces the number density of sources to no = 
0.0087 Mpc^^(30 kpc/r)^, which is similar to the density of galaxies brighter than 0.12L*. 
In other words, if every L > O.IL* galaxy goes through a starburst phase, and remnants of 
the cool outflows persist into the post-starburst phase, then relict winds could account for 
the entire path density of strong Mg II systems at low redshift. Larger cross sections clearly 
raise the minimum mass of galaxies that would need to contribute relic outflows to the path 
density. 

This hypothesis needs to be tested further against two additional constraints. First, 
the local path density of OVI systems, which may trace the shocks from these winds, also 
suggests dwarf galaxies make a large contribution to the O VI cross section (Danforth & 
Shull 2005; Stocke etal. 2005). However, to date, direct absorber - galaxy associations 
require cross sections (Tr ~ 200 kpc (Adelberger etal. 2005; Stocke etal. 2005; Tumlinson & 
Fang 2005), which is plausible but not yet required for the cool outflows. Second, and more 
importantly, the velocity widths of low-ionization, intervening systems need to be compared 
to the distribution measured for outflows, which has a mean here of 350 km s^^ FWHM. 
Songaila et al. (2005b) suggest that the linewidths and ionization state of absorbers can 
differentiate winds from other types of halo gas. 



6. Conclusions 

This work maps the outflow properties across 18 ultraluminous infrared-selected galax- 
ies, providing measurements of Na I column density, Na I absorption-line velocities, and Ha 
emission-line velocities. These starburst galaxies are undergoing major mergers, as evidenced 
by the presence of multiple nuclei and prominent tidal features. The cool outflows previously 
identifled in absorption against the stellar continuum (Heckman et al. 2000; Rupke et al. 2002; 
Martin 2005; Rupke et al. 2005a, 2005b) are shown to cover a large portion of the galaxy. 
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exceeding 15 kpc along the slit in a few cases. Faint, spatially-extended Ha emission, often 
diffuse emission or narrow spikes in the spectra rather than well-defined Doppler ellipses as 
found in IRAS0315+42, is detected in about half of the ULIGs; but its presence shows no 
obvious correlation with the strength of the Na I absorption. The gas velocity in the blue 
wing of the Ha line profile is comparable to, and may exceed (Rupke etal. 2005b), that of 
the most blue-shifted absorbing gas. It follows that the outflow probably does contain warm 
Ha -emitting gas; but the emission line profile is largely shaped by higher density gas. 

The absorption profiles extend from the systemic velocity up to (blueshifted) veloci- 
ties comparable to the halo escape velocity, suggesting acceleration of cool, interstellar gas. 
Across many ULIGs, IRAS1056-I-24 being one example, the outflow velocity is constant over 
the entire galaxy while the Ha velocity shows a marked rotational gradient, shifting from 
negative to positive velocities across the nucleus. These results are consistent with an outflow 
diverging from a central starburst, because the rotational velocity of a diverging wind slows 
down with increasing height to conserve angular momentum. Surprisingly, however, across 
several ULIGs, particularly IRAS1720-00 and IRAS0015+54, the cool outflows present a ve- 
locity shear as large as the projected galactic and orbital rotation combined. These marked 
kinematics suggest a low scale height and large launch region for at least these outflows. It 
will be interesting to see whether numerical simulations of winds from galaxy mergers can 
reproduce such results by including spatially extended star formation or merger-driven out- 
flows. If not, coeval winds from both nuclei might offer the only explanation as overlapping 
disks do not work in general. 

The Na I AA5890, 5896 lines are blended in ULIG spectra, due to the presence of mul- 
tiple absorption components along the sightline. Where this absorption trough was highly 
structured, the doublet ratio yielded optical depths ~ 2 to 10. Less structured troughs 
allowed much larger optical depths, which is problematic due to the nonlinear relation be- 
tween the doublet ratio and the column density. This degeneracy was resolved by fitting 
the line optical depth, rather than the doublet ratio. The fractional area of the beam cov- 
ered by absorbing gas was estimated from the residual line intensity before calculating the 
optical depth profiles from the intensity profiles. This approach produced column densities, 
iV(Na I), between 1.8 x 10^^ cm~^ and 2.4 x 10^^ cm~^. A few systems common to the Rupke 
et al. (2005a) sample, and our Na I column density measurements agree to within a factor 
of two. 

The mass of outfiowing gas, estimated by integrating column densities over the area of 
the outfiow, is of order 10^ Mq, or a few percent of the gas mass, in most ULIGs. These 
values are lower limits in several systems where the projected outfiow appears to be larger 
than the background continuum source. The reddening increases with the strength of the 
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Na I absorption as was found previously in LIGs (Heckman et al. 2000) but exceeds the 
amount expected for the estimated outflow columns and a Galactic dust-to-gas ratio. The 
applied corrections for Na ionization and depletion could be low, indicating an even higher 
gas column, but this is not necessary as some of the dust may not make it into the wind, 
just as much of the interstellar gas is left behind. Previous empirical results suggested mass 
outflow rates in the hot and warm phases were comparable to the SFR in starburst galaxies 
(Martin 1999), although this is not well-constrained in galaxies as luminous as ULIGs. Based 
on a simple model for a freely-expanding outflow, the measured gas columns and velocities 
indicate mass loss rates in the cool gas of order 10% the SFR, in rough agreement with Rupke 
etal. (2005b). This efficiency would be higher in dwarf starbursts than ULIGs if ionization 
and depletion corrections are actually larger in the former as suggested. Such uncertainties 
are currently too large to allow a definitive discrimination between momentum-driven and 
energy-driven outfiows in the cool gas. 

The fate of the cool outfiows remains model dependent for now. Merger-driven winds 
appear unlikely to eject a lot of material from the halo potential, and their viability is not 
yet certain. It is clear that momentum- driven winds would easily transport the gas detected 
in Na I absorption to radii ~ 30 kpc, and much further for some choices of parameters. The 
momentum carried by the cool winds is roughly 10% of that supplied by supernovae and 
stellar winds (or the radiation field); and the kinetic energy carried by the cool gas is a few 
percent of the supernova energy (for a burst age of 10 Myr). These measured efficiencies 
provide a guide for feedback recipes in galaxy formation models, which should distinguish 
the cool and hot gas. They appear consistent with estimates for LIGs and ULIGs presented 
by Rupke et al. (2005b) and theoretical expectations for supernova-driven winds. 

The interpretation of the data presented here could be improved by new numerical sim- 
ulations. Numericists are beginning to address the temporal evolution of outfiows (Hopkins 
etal. 2005), spatially extended star formation (Fragile etal. 2004), and merger-driven winds 
(Cox etal. 2004; Cox etal. 2005). It will be interesting to learn what rotational properties 
these models generate for the entrained gas component, and whether clouds dense enough 
to produce Na I absorption survive in the wind. More detailed simulations of merger-driven 
winds are needed to determine their viability, particularly quantitative predictions for the 
efficiency with which orbital kinetic energy is transferred to the bulk motion of the outfiow. 
More observations of thermal X-ray emission from ULIGs will distinguish whether a hot 
wind plays an important role in the acceleration of the cool gas, as the hot phase is not 
required in a purely radiative model of acceleration. 

The measurements presented here suggest the cross sections and column densities for 
the cool component of winds are significant. Ultraluminous infrared galaxies are rare in the 
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local universe, but such extremely luminous starbursts were much more common in the past. 
If most galaxies brighter than 0.1 L* have undergone such extreme star formation activity, 
the relic winds should be detectable as intervening absorption-line systems. 

The author thanks TJ Cox, Tim Heckman, Norm Murray, Todd Thompson, Eliot 
Quataert, and Sylvain Veilleux for stimulating discussions. The R band images shown in 
this paper were kindly provided by Lee Armus and Tom Murphy. The author thanks Dr. 
Colleen Schwartz for flux calibrating the spectra. This paper was written largely while the 
author was a visitor at the KITP. This work was supported by grants from the David and 
Lucile Packard Foundation and the Alfred P. Sloan Foudation. 

A SIMPLE TOY MODEL FOR THE ROTATION OF CONICAL WINDS 

When clouds launched from a galactic disk follow ballistic trajectories through the halo, 
they rotate more slowly than the gas disk (Bregman 1980; Collins etal. 2002). A starburst- 
driven wind differs from these fountain models in that the starburst outflow is launched from 
a much smaller region of the disk. The expansion of the superbubble quickly stalls in the 
dense gas disk, and the shell develops a bipolar shape (De Young & Heckman 1994; Strickland 
& Stevens 2000; MacLow & Ferrara 1999; Mori etal. 2002, Fragile etal. 2004). After the 
elongated supershell erupts from the galaxy, the outflow is mostly radial (Chevalier & Clegg 
1985). Rotation of the outflow has been ignored in these calculations largely because it is 
expected to be small. 

To provide a baseline model for the rotational properties of galactic outflows, as mea- 
sured by absorption-line position-velocity diagrams, this Appendix considers the rotation of 
galactic winds in two limits: (1) constant angular momentum, and (2) constant angular veloc- 
ity. Angular momentum conservation provides the simplest model. Outflows with constant 
angular velocity provide some insight into the rotational properties of magnetically-driven 
winds. Only when such torques act on the outflow, will significant rotation be observed in 
projection along the major axis of the disk. When the velocities are projected against the 
minor axis, a velocity jump can be seen across the nucleus, in the absence of any torques. 
The argument goes as follows. 



A. Al. OUTFLOWS CONSERVING ANGULAR MOMENTUM 

Consider the trajectory of clouds launched from a central starburst in a disk, which 
rotates at speed Vc- Figure 11 sketches the geometry. The starburst region is described by 
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a squat cylinder of radius Rq. Disk gas entrained in the wind at Rq flows freely outward 
at speed V^, and this outflow diverges radially along a conical surface of opening angle 
6u,- The thickness of the starburst region is Zq = Ro/tan 6w As a ring of wind material 
moves outward along the cone, the circular velocity of the gas declines to conserve angular 
momentum as V{R) = VcRo/R, where R is the radius in cylindrical polar coordinates. 

Now consider what an observer viewing this simple toy model would see. The location of 
the observer in Figure 11 deflnes the inclination of the galaxy, i, which is the angle between 
the sightlines and the polar axis of the galaxy. A slit aligned with the major axis of this 
galaxy samples the model along a series of parallel sightlines that intersect the x-axis (see 
Figure 11). The intersection of these sightlines with the surface of the conical outflow deflne 
the locus where the velocity fleld is sampled. In each plane, x = constant, the intersection 
of the sightline and the hyperbolic conic section is easily computed. At these points, the 
projection of the circular velocity along the sightline is 

Vcproj = V{R)cos{(f))sin{L), (Al) 

and the projection of the outflow velocity is 

Vw^proj = VwSin{6w)sin{(l))sin{L) + VwCos(6yj)cos{L). (A2) 

For a wind half-opening angle of ^^ = 60° , the resulting position-velocity diagram for 
projection along the major axis is shown in Figure 12a. The rotation speed of the disk 
changes sign on the approaching and receding sides of the major axis, but the projected 
speed of the radial wind is symmetric to either side of the minor axis. Along the major axis 
of a galaxy, the outflow velocities measured in absorption will not present much of a velocity 
gradient compared to the rotation curve of the disk provided V^ > Vc- 

For sightlines intersecting the minor axis, the rotational velocity is perpendicular to 
the sightline and therefore undetectable. Projection of the radial outflow gives observed 
velocities 

yw.pTO] = yw{,sin{Bu^s%n{ji) ^ cos{dy^cos{i)\y^ (A3) 

= Vw{^—sin{Qw^sin{i) + cos(Qw)cos{i)\y < 0, (A4) 

along the minor axis (y-axis in Figure 11). The P-V diagram along the minor axis is not 
symmetric about the nucleus. On the near side of the minor axis, where the cone is tipped 
toward us, the apparent velocity may approach \4). In contrast, on the far side, the radial 
velocity vector makes a larger angle with the sightline, so the projected speed is lower. 
Figure 12b illustrates the jump in the projected wind speed along the minor axis for several 
viewing angles. Notice that the wind can appear redshifted against the far side of the disk 
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when the viewing angle i is close to the opening angle of the cone - i.e. 90° —6w<i< Oy,. In 
summary, although the wind velocity varies little along the major axis, it can jump abruptly 
along the minor axis near the nucleus. The P-V diagram would not show a smooth velocity 
gradient along the entire minor axis. 



B. A2. CONSTANT ANGULAR VELOCITY WIND 

This physically-motivated, but purely geometrical, model provides some intuitive feel for 
how projection affects the observed wind velocities. Combined with observations, it suggests 
some outflows acquire additional angular momentum as they leave the disk. Magnetic field 
lines anchored in the disk but rising into the outflow cone would act somewhat like rigid wires 
because the gas is forced to corotate with the field lines out to the Alfven point (Pino & 
Tanco 1999). Magnetic fields might therefore spin up the outflow at the expense of the disk 
angular momentum. To illustrate the potential impact on the observed velocity gradient, 
the toy model was given a constant angular velocity. The circular velocity at every point was 
set to V,j) = VcR/Ro- The wind velocity is again radial and was taken to be approximately 
constant at V^, ignoring magnetic slingshot effects. The blue lines in Figure 12a illustrate 
the projected P-V diagram along the major axis. The projected velocity gradient is much 
steeper than the constant-L model. For a wind speed twice that of the disk rotation, the 
outflow appears blueshifted, toward the receding side of the major axis, out to two starburst 
radii. This example demonstrates that outflows with constant angular velocity can produce 
blueshifted absorption with large velocity gradients. The presence of halo magnetic fields is 
supported by observations of some starbursts, which show extended synchrotron emission 
(e.g. Thompson etal. 2006; Israel & deBruyn 1988; Dahlem etal. 1995). 



-32- 



REFERENCES 

Adelberger, K. L. etal. 2005, ApJ, 629, 636 

Arav, N. et al. 1999b, ApJ, 524, 566 

Arav, N. etal. 1999a, ApJ, 516, 27 

Annus, L., Heckman, T. M., & Miley, G. K. 1989, ApJ, 347, 727 

Barlow, T. A. & Sargent, W. L. W. 1997, AJ113, 136 

Bendo, G. J. & Barnes, J. 2000, MNRAS, 316, 315 

Blanton, M. etal. 2003, ApJ, 592, 819 

Bregman, J. N. 1980, ApJ, 236, 577 

Chevalier, R., & Clegg, A. W. 1985, Nature, 317, 44 

Colina, L. etal. 2004, ApJ, 602, 81 

Collins, J. A., Benjamin, R. A., & Rand, R. J. 2002, ApJ, 578, 98 

Cox, T. J. etal. 2005, astro-ph-0504156 

Cox, T. J. etal. 2004, ApJ, 607, L87 

Dahlem, M., Lisenfeld, U., & Golla, G. 1995, ApJ, 444, 119 

Danforth, C. W. & Sliull, J. M. 2005, ApJ, 624, 555 

DeYoung, D. S., & Heckman, T. M. 1994, ApJ, 431, 598 

Downes, D., Solomon, R M., & Radford, S. J. E. 1993, ApJ, 414, L13 

Ferlet, R., Vidal-Madjar, A., & Cry, C. 1985, ApJ, 298, 838 

Fragile, R C., Murray, S. D., & Line, D. N. C, astro-ph/0408438 

Franx, M. etal. 1997, ApJ, 486, 75 

Fujita, A. et al. 2006, in prep 

Hamann, F., Barlow, T. A., Junkkarinen, V., & Burbidge, E. M. 1997, ApJ, 478, 80 

Heckman, T., Lehnert, M. D., Strickland D. K., & Lee, A. 2000, ApJS, 129, 493 

Heckman, T. M., Armus, & Miley 1990 

Heckman, T. M. & Leitherer, C. 1997, A J, 114, 69 

Genzel, R. etal. 2001, ApJ, 563, 527 

Hibbard, J. E. etal. 2000, AJ, 119, 1130 

Hopkins, R F. 2005, ApJ, 630, 705 



-33- 



Hunter, D. A., Thronson, H. A., Jr. 1996, ApJ, 461, 202 

Israel, F. P. & de Bruyn, A. G. 1988, A&A, 198, 109 

Jenkins, E. O. 1986, ApJ, 304, 739 

Jog, C. J. & Solomon, P. M. 1992, ApJ, 387, 152 

Johnson, K. etal. 2000, AJ, 120, 1273 

Kriss, G. 1994, in ASP Conf. Ser. 61, Astronomical Data Analysis Software and Systems 
III, ed. D. R. Crabtree, R. J. Hanisch, & J. Barnes (San Francisco: ASP), 437 

Kunth, D. etal. 1998, A&A, 334, 11 

Lehnert, M. D., & Heckman, T. M. 1996, ApJ, 462, 65 

Leitherer, C. etal. 1999, ApJS, 123, 3 

MacLow, M.-M., & Ferrara, A. 1999, ApJ, 513, 142 

Martin, C. L. 2005, ApJ, 621, 227 (Paper I) 

Martin, C. L. 1998, ApJ, 506, 222 

Massey, P. & Gronwall, C. 1990, ApJ, 358, 344 

McDowell, J. etal. 2003, ApJ, 591, 154 

Mihos, J. C. 1999, Ap&SS, 266, 195 

Miller, J. S., & Mathews, W. G. 1972, ApJ, 172, 593 

Morton, D. 1991, ApJS, 77, 119 

Mori, M. Ferrara, A., & Madau, P. 2002, ApJ, 571, 40 

Murphy, T. W., Jr., Soifer, B. T., Matthews, K., & Armus, L. 2001, ApJ, 559, 201 

Murphy, T. W., Jr. etal. 1996, AJ, 111, 1025 

Murray, N. etal. 2006, ApJ, submitted 

Murray, N., Quataert, E., & Thompson, T. A. 2005, ApJ, 618, 569 

Pettini, M. etal. 2002 ApJ, 569, 742 

Phillips, A. C. 1993, AJ, 105, 486 

Pino, E. & Tanco, G. 1999, ApJ, 518, 129 

Prochter, G. E., Prochaska, J. X., & Buries, S. M. 2004, astro-ph/0411776 

Puget, J.-L. etal. 1996, A&A, 308, L5 

Rand, R. J. 2000, ApJ, 537, L13 



-34- 



Rupke, D. S., Veilleux, S., & Sanders, D. B. 2002, ApJ, 570, 588 

Rupke, D. S., Veilleux, S., & Sanders, D. B. 2005a, ApJS, 160, 87 

Rupke, D. S., Veilleux, S., & Sanders, D. B. 2005b, ApJS, 160, 115 

Rupke, D. S., Veilleux, S., & Sanders, D. B. 2005c, ApJ, 631, L37 

Keck High-Resolution Spectroscopy of Outflows in Infrared-luminous Galaxies 

Rupke, D., Veilleux, S., & Sanders, D. B. 2005d, ApJ, 632, 751 

Sanders, D. B. 1988, ApJ, 325, 74 

Savage, B. D. & Sembach, K. R. 1996, ARAA, 34, 279 

Savage, B. D. etal. 1977, ApJ, 216, 291 

Schwartz, C. S. & Martin, C. L. 2004, ApJ, 610, 20 

Schwartz, C S. 2005, PhD Thesis, Univ. of California, Santa Barbara 

Scoville, N. Z. etal. , 2000, ApJ, 119, 991 

Scoville, N. Z., Sargent, A. I., Sanders, D. B., & Soifer, B. T. 1991, ApJ, 366, 5 

Shapley, A. E. et al. 2003, ApJ, 588, 65 

Sheinis, A. I. etal. 2002, PASP114, 851 

Soifer, B. T. etal. 1986, ApJ, 303, L41 

Soifer, B. T. etal. 1987, ApJ, 320, 238 

Solomon, P. M., Downes, D., Radford, S.J.E, & Barrett, J. W. 1997, ApJ, 478, 144 

Songaila, A. 2005a, astro-ph/0507649 

Songaila, A. 2005b, astro-ph/0509821 

Songaila, A. 2001, ApJ, 561, 153 

Spitzer, L. 1968, Diffuse Matter in Space (New York: Interscience) 

Strauss, M. A. etal. 1992, ApJS, 83, 29 

Steidel etal. 2002, ApJ, 570, 526 

Steidel, C. S. & Sargent, W. W. 1992, ApJS, 80, 1 

Stocke, J. T. etal. 2005, astro-ph/0509822 

Strickland, D. K. & Stevens, I. R. 2000, MNRAS, 314, 511 

Strickland, D. K. etal. 2004a, ApJ, 151, 193 

Strickland, D. K. etal. 2004b, ApJ, 606, 829 



-35- 



Surace, J. A. & Sanders, D. B. 2000, AJ, 120, 604 

Tacconi, L. J. etal. 2002, ApJ, 580, 73 

Thompson, T., Quartaert, E., Waxman, E., Murray, N., & Martin, C L. etal. 2006, astro-ph 

Tumlinson, J. & Fang, T. 2005, ApJ, 623, L97 

Tullmann, R. et al. 2000, A&A, 364, L36 

Vazquez, G. A. etal. 2004, ApJ, 600, 162 

Veilleux, S. etal. 1995, ApJS, 98, 171 

Veilleux, S., Cecil, G., & Bland-Hawthorn, J. 2005, ARA&A, 43, 769 

Wakker, B. R & Mathis, J. S. 2000, ApJ, 544, L107 

Weaver, B. etal. 1977, ApJ, 218, 377 

WeiB, A. et al. A&A, 365, 571 



This preprint was prepared with the AAS I^TJtjX macros v5.2. 



-36- 



Table 1: Properties of Extended Na I Absorption 
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2'.'6 


6.2 


4.7 X 10^ 


5.6 X 


1055 


11 


1.2 X 10^ 


IRAS03521+0028 


76.5 


1'.'9 


5.0 












IRAS08030+5243 





>5'.'2 


>8.2 


>1.6 X 10^ 






3 


7.7 X 10^ 


IRAS 10494+4424 


25.0 


>7^0 


>12 


>5.9 X 10^ 


1.4 X 


1055 


>3 


1.2 X 10^ 


IRAS10565+2448 


109.0 


>13'.'4^ 


>11 


>1.1 X 10^ 


4.6 X 


1055 


>8 


2.3 X 10^ 


IRAS11506+1331 


80.3 


5'.'9^ 


13 


1.74 X 10^ 


4.3 X 


1054 


0.72 


3.2 X 10^ 


IRAS15245+1019 


127.8 


>12'.'4 


>17.8 


>6.7 X 10^ 


2.8 X 


1055 


>6 


4.4 X 10^ 


IRAS16090-0139 


107.6 


3'.'2 


7.6 


8.0 X 10^ 


7.0 X 


1055 


8 


1.8 X 10^ 


IRAS16487+5447 


66.2 


5'.'7 


11 












IRAS17208-0014 


166.7 


>8'.'6 


>7.3 


7.0 X 10^ 


3.8 X 


1055 


7 


1.2 X 10^ 


IRAS18368+3549 


-31.0 


4'.'8 


10 


1.4 X 10« 


1.2 X 


1056 


5 


1.9 X 10^ 


IRAS19297-0406 


149.5 


2'.'9 


4.8 


4.0 X 10^ 


5.8 X 


1055 


5 


5.3 X 10^ 


IRAS19458+0944 


118.1 
















IRAS20087-0308 


85.5 


>6'.'3 


>12 


>9.5 X 10^ 


1.1 X 


1056 


>14 


3.8 X 10^ 


IRAS23365+3604 


-19.5 


4'.'7 


5.9 


4.4 X 10^ 


5.6 X 


1055 


8 


1.1 X 10^ 





















Note. — Col 1 - Galaxy name; see Table 1 of Paper I for redshift and star formation rate; Col 2 - 
Slit position angle (east of north); Col 3 - Angular length of slit where Na I absorption detected. Lower 
limits indicate detections extend as far as continuum detection. Col 4 - Projected spatial extent of the Na I 
absorbing gas, scales ash^Q] Col 5 - Mass of gas detected in absorption. The corrections for dust depletion 
and ionization correction are as given in Eqn. 3, and the Na I columns from Table 3 column 6 are adopted. 
The area for each aperture is assumed to be an annulus centered on the nucleus, and the masses of all 
apertures are summed. Col 6 - Kinetic energy estimated from the summed apertures KE — ^ MiV^^, using 
the same masses as in Col. 6 and the velocities from Col. 9 of Table 3. Col 7 - Mass loss rate estimated 
from a free-streaming wind model normalized to the column density and velocity in the nuclear aperature, 
as described by Eqn. 7. Col 8 - Mass in wind estimated from the model described by Eqn. 8. The launch 
radius is taken from Col. 3 Table 2 , and the maximum extent is taken from Col. 4 above. 
"Detected Na I toward two nuclei, and the entire extent of the absorption is listed. 
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Table 2: Emission Properties 



Object 


Scale 


Rco 


Ha Extent 


VvHa (km 


VvNai (km 


c(H/3) 


Merge 




(kpc/" ) 


(pc) 


(kpc) 


s^i kpc-i) 


s^i kpc-1) 




State 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


IRAS00153+5454 


2.038 




25.5 (blob) 


7 (17, 12 & 7) 


16 + 4 


1.41 + 0.82 


doubh 


IRAS00188-0856 


2.295 


297 


10 (halo) 





7 + 16 


2.42 + 2.42 


doubh 


IRAS00262+4251 


1.799 


314 


19.8 (halo) 


15 


7+19 


0.99 + 0.37 


tails 


IRAS03158+4227 


2.384 


405 


31.6 (spike) 


16 


40 + 14 


1.72 + 3.05 


single 


IRAS03521+0028 


2.641 


575 








2.81 + 2.81 


irreg 


IRAS08030+5243 


1.570 


345 


19.9 (blobs) 


60 


8 + 3 


1.33 + 0.45 


single 


IRAS 10494+4424 


1.717 


379 


17.6 (spike) 


2 


24 + 5 


1.79 + 3.83 


tails 


IRAS10565+2448 


0.849 


317 


6.4 {r.h) 


13 (11, 19) 


9 + 7 


1.45 + 0.53 


doubh 


IRAS11506+1331 


2.276 


361 


27.7 (spike) 


3 (1, 14 & 7) 


10 + 5 


1.23 + 1.30 


doubh 


IRAS15245+1019 


1.436 




18.0 (spike) 


9 (10, 6 & 21) 


1 + 1 


1.47 + 0.69 


doubh 


IRAS16090-0139 


2.371 


459 


20 (spike) 


8 


19 + 19 


0.99 + 2.93 


single 


IRAS16487+5447 


1.903 




11.6 (diffuse) 






0.41 + 0.41 


doubh 


IRAS17208-0014 


0.846 


301 


7.9 


47 (38, 65 & 20) 


42 + 2 


1.53 + 2.50 


doubh 


IRAS18368+3549 


2.104 


396 


21.2 (halo) 


35 


23 + 11 


2.6 + 2.6 


tails 


IRAS 19297-0406 


1.607 


385 


15 


24 


2 + 9 


1.78 + 3.83 


single 


IRAS 19458+0944 


1.844 


465 








1.79 + 1.79 


tails 


IRAS20087-0308 


1.936 


466 


6.5 


26 (15, 25 & 5) 


19 + 8 


4.0 + 4.0 


doubh 


IRAS23365+3604 


1.239 


372 


15.4 (blobs) 


34 (0, 34) 


6 + 29 


1.20 + 0.50 


tails 



















Note. — (1) Object; (2) Angular scale at redshift given in Paper I (h = 0.7, ilo = 0.3, ^a = 0.7); (3) 
Radius of CO emission from Solomon et al. (1997); (4) Spatial extent of Ha emission along slit (PA given in 
col. 2 of Table 1, and slit center shown in Figure 1); and description of Ha morphology based on the upper 
right panels of Figure 1. (5) Total Ha velocity gradient across the ULIG. For double systems, the projected 
orbital and rotational components are listed, respectively, in parentheses. (6) Total Na I velocity gradient 
across the ULIG. (7) Logarithmic extinction at H/3 dervied from the measured Ha /H/3 flux ratio, an intrinsic 
ratio of 2.86, and a standard interstellar extinction curve (Miller & Mathews 1972). The difference between 
this value and that derived from the F(H/3 )/F(H7 ) flux ratio was used to estimate the systematic error in 
the reddening. (8) Merger morphology based on Murphy et al. (2001) scheme and images from Murphy ct 
al. (1996). 



Table 3. Na I Spectral Line Fitting 



IRAS 


Ap 


Offset 


Cf 


r 


N(Nal) 


DR 


N(Nal) 


V 


AV 


min N(Nal) 


max N(Nal) 


V 


AV 


Vhc 






(") 




(A5896) 


(cm 


-^) 




(cm 


-^) 


(kms-1) 


(km s-i) 


(cm-2) 


(cm-2) 


(km s-i) 


(km s-i) 


(km s-i) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


(15) 


0015+54 


1 


2.39 




















6.3 X 


10" 


3.3 X 10" 


-156± 191 


449+ 391 


-150 


0015+54 


2 


0.70 


0.30 


6.88 


7.0 X 


10« 


1.137 


4.5 X 


1013 


-364 ± 19 


555 ± 54 


4.3 X 


1013 


5.5 X 10" 


-372± 21 


533+ 55 


-143 


0015+54 


3* 


-0.98 


0.35 


4.79 


4.9 X 


IQl-^ 


1.167 


2.6 X 


1013 


-276 ± 34 


497 ± 89 


2.1 X 


1013 


4.2 X 10" 


-285± 34 


462+ 85 


-65 


0015+54 


4 


-4.2 


0.58 


4.31 


4.4 X 


1013 


1.177 


4.1 X 


1013 


-183 ± 20 


452 ± 48 


5.1 X 


1013 


7.4 X 10" 


-199+ 22 


442+ 47 


28 


0015+54 


5 


-5.9 


0.40 


3.12 


3.2 X 


10" 


1.219 


3.3 X 


1013 


-179 ± 21 


450 ± 51 


1.3 X 


1013 


3.6 X 10" 


-158+ 52 


380+ 110 


-44 


0018-08 


1 


6.99 




















9.1 X 


109 


9.1 X 109 


315± 302 


585+ 1214 




0018-08 


2 


1.64 


0.24 


5.55 


5.7 X 


1013 


1.153 


1.8 X 


1013 


-211 ± 141 


421 ± 130 


7.5 X 


10" 


2.8 X 10" 


-214± 175 


414+ 140 


-16 


0018-08 


3* 





0.45 


6.53 


6.7 X 


1013 


1.141 


5.8 X 


1013 


-356 ± 13 


574 ± 19 


1.2 X 


1014 


8.0 X 10" 


-358+ 13 


573+ 18 


-33 


0018-08 


4 


-1.72 


0.17 


7.33 


7.5 X 


1013 


1.133 


2.7 X 


1013 


-355 ± 85 


609 ± 104 


1.1 X 


1013 


3.3 X 10" 


-355+ 100 


611+ 109 


-17 


0026+42 


1* 


0.48 


0.14 


4.24 


2.4 X 


IQi-i 


1.179 


1.3 X 


1013 


-211 ± 18 


387 ± 35 


2.1 X 


1013 


2.4 X 10" 


-218+ 18 


384+ 36 


22 


0026+42 


2 


-1.19 


0.12 


4.84 


1.7 X 


1014 


1.163 


9.3 X 


10" 


-284 ± 21 


357 ± 37 


1.2 X 


1013 


1.7 X 10" 


-293+ 21 


401+ 49 


-53 


0026+42 


3 


-2.88 




















9.2 X 


IQii 


4.2 X 10i3 


-263+ 116 


281+ 124 


-39 


0315+42 


1 


1.68 




















3.7 X 


10" 


6.1 X 10" 


-708+ 128 


670+ 221 


49 


0315+42 


2* 





0.41 


7.00 


7.2 X 


1013 


1.135 


5.7 X 


1013 


-526 ± 14 


461 ± 49 


9.6 X 


10" 


1.0 X 10i3 


-465+ 18 


514+ 45 


37 


0315+42 


3 


-1.68 


0.28 


7.16 


7.3 X 


1013 


1.135 


3.9 X 


1013 


-322 ± 28 


381 ± 61 


5.9 X 


10" 


7.3 X 10" 


-306+ 51 


426+ 66 


-27 


0803+52 


1 


1.94 


0.69 


3.84 


3.9 X 


1013 


1.194 


1.8 X 


1013 


238 ± 4 


160 ± 6 


1.0 X 


1013 


4.5 X 10" 


233+ 9 


164+ 16 


-252 


0803+52 


2* 





0.76 


4.05 


4.1 X 


1013 


1.187 


3.1 X 


1013 


258 ± 3 


150 ± 4 


1.8 X 


1013 


5.7 X 10i3 


259+ 2 


147+ 4 





0803+52 


3 


-1.94 


0.48 


4.18 


4.3 X 


1013 


1.184 


5.5 X 


10" 


276 ± 13 


187 ± 17 


4.9 X 


10" 


1.0 X 1013 


281+ 13 


184+ 17 


115 


1049+44 


1 


3.00 


0.57 


7.22 


7.9 X 


1013 


1.132 


9.0 X 


1013 


-339 ± 40 


539 ± 66 


7.6 X 


1013 


1.0 X IQis 


-347+ 21 


532+ 59 


-3 


1049+44 


2 


1.37 


0.57 


7.83 


8.0 X 


1013 


1.131 


7.6 X 


1013 


-227 ± 43 


475 ± 36 


4.1 X 


1013 


9.3 X 10" 


-228+ 43 


472+ 35 


24 


1049+44 


3* 


-0.27 


0.65 


3.70 


3.8 X 


1013 


1.198 


1.1 X 


1013 


-300 ± 39 


273 ± 38 


1.1 X 


1013 


2.5 X 10" 


-297+ 37 


271+ 36 


-65 


1049+44 


4 


-2.35 


0.40 


2.80 


2.9 X 


1013 


1.241 


1.1 X 


1013 


-114 ± 41 


374 ± 


8.4 X 


10" 


2.6 X 10" 


-122+ 41 


237+ 56 


-17 


1056+24 


1 


2.98 


0.80 


1.97 


2.2 X 


1013 


1.310 


1.7 X 


1013 


-427 ± 35 


329 ± 57 


2.6 X 


1013 


1.2 X lOis 


-446+ 33 


332+ 56 


52 


1056+24 


2* 


0.45 


0.78 


6.89 


7.0 X 


1013 


1.142 


5.2 X 


1013 


-466 ± 8 


249 ± 10 


1.1 X 


1014 


1.6 X 10i5 


-460+ 7 


294+ 9 


57 


1056+24 


3 


-1.04 


0.86 


4.20 


4.3 X 


1013 


1.190 


3.2 X 


1013 


-392 ± 190 


403 ± 207 


3.4 X 


1013 


1.5 X lOis 


-394+ 36 


437+ 42 


-41 


1056+24 


4 


-4.17 


0.80 


1.78 


1.8 X 


1013 


1.353 


8.5 X 


109 






2.1 X 


10" 


7.7 X 10" 


-467+ 274 


487+ 210 


-64 


1056+24 


5 


-8.34 


0.47 


1.89 


2.0 X 


1013 


1.337 


9.0 X 


10" 


-358 ± 67 


467 ± 126 


9.7 X 


10" 


6.1 X 10" 


-356+ 58 


388+ 170 


-70 


1056+24 


lb 


2.98 


0.80 


3.75 


3.7 X 


1013 


1.208 


2.9 X 


1013 


-210 ± 6 


138 ± 13 


4.2 X 


1013 


6.5 X 10" 


-211+ 5 


145+ 11 




1056+24 


2b* 


0.45 


0.78 


2.22 


2.3 X 


1013 


1.298 


2.8 X 


1013 


-186 ± 7 


208 ± 12 


7.3 X 


1013 


7.3 X 10" 


-194+ 6 


242+ 9 




1056+24 


3b 


-1.04 


0.86 


3.39 


3.5 X 


1013 


1.216 


3.6 X 


1013 


-151 ± 3 


169 ± 8 


5.7 X 


1013 


8.2 X 10" 


-150+ 2 


191+ 11 




1056+24 


4b 


-4.17 


0.80 


2.44 


2.5 X 


1013 


1.284 


2.3 X 


1013 


-181 ± 5 


125 ± 9 


4.6 X 


1013 


6.7 X 10" 


-181+ 4 


153+ 12 




1056+24 


5b 


-8.34 


0.47 


2.91 


3.0 X 


1013 


1.246 


1.2 X 


1013 


-113 ± 7 


148 ± 7 


1.1 X 


1013 


3.0 X 10" 


-108+ 4 


163+ 14 




1150+13 


1 


5.51 




















2.3 X 


10" 


9.4 X 10i3 


66+ 77 


285+ 103 


22 



00 
00 























Table 3 — Continued 














IRAS 


Ap 


Offset 


Cf 


r 


N(Nal) 


DR 


N(Nal) 


V 


AV 


min N(Nal) 


max N(Nal) 


V 


AV 


Vhc 








(") 




{A5896) 


(cm 


-') 




(cm 


-') 


(km s-1) 


(km s-i) 


(cm- 2) 


(cm-2) 


(km s-i) 


(km s-i) 


(km s-i) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 




(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


(15) 




1150+13 


2 


4.05 


0.21 


4.22 


4.3 X 


lOi^ 


1.179 


1.0 


X 


1013 


63 ± 30 


245 ± 34 


9.0 X 


1012 


1.9 X 10" 


56± 32 


247+ 36 


-89 




1150+13 


3 


1.53 


0.22 


2.12 


2.2 X 


10i« 


1.293 


1.7 


X 


1012 


139 ± 74 


221 ± 


5.4 X 


1012 


2.6 X 10" 


-43± 55 


433+ 94 


73 




1150+13 


4* 


-0.15 


0.43 


3.10 


3.2 X 


1013 


1.219 


2.2 


X 


10" 


-89 ± 3.5 


221 ± 4 


2.7 X 


10" 


5.2 X 10" 


-90± 3 


224+ 4 


8 




1150+13 


5 


-1.84 


0.33 


2.72 


2.8 X 


1013 


1.243 


1.0 


X 


10" 


-104 ± 14 


157 ± 19 


1.1 X 


10" 


3.0 X 10" 


-103+ 14 


160+ 18 


-35 




1524+10 


1 


3.61 


0.59 


2.74 


2.8 X 


1013 


1.247 


1.7 


X 


10" 


-200 ± 3 


138 ± 7 


1.4 X 


10" 


2.5 X 10" 


-201+ 2 


127+ 3 


-144 




1524+10 


2 


1.52 


0.61 


6.67 


6.8 X 


10" 


1.142 


6.7 


X 


10" 


-215 ± 5 


312 ± 12 


3.2 X 


10" 


5.1 X 10" 


-211+ 5 


297+ 7 


-77 




1524+10 


3* 


-0.72 


0.78 


7.11 


7.3 X 


10" 


1.138 


9.9 


X 


10" 


-245 ± 1 


310 ±4 


8.0 X 


10" 


2.6 X 10" 


-247+ 3 


307+ 5 


-75 




1524+10 


4 


-2.35 


0.84 


5.56 


5.6 X 


10" 


1.157 


7.3 


X 


1013 


-195 ± 4 


260 ± 13 


6.5 X 


10" 


4.1 X 10" 


-194+ 2 


217+ 5 


-42 




1524+10 


5 


-4.14 


0.75 


3.85 


3.9 X 


10" 


1.194 


4.0 


X 


10" 


-186 ± 7 


218 ± 16 


3.0 X 


10" 


7.5 X 10" 


-188+ 5 


183+ 9 


-43 




1609-01 


1 


3.37 


0.58 


2.75 


2.8 X 


10" 


1.239 


1.3 


X 


10" 


-314 ± 25 


132 ± 32 


3.6 X 


1012 


6.9 X 1012 


-328+ 22 


116+ 31 


-4 




1609-01 


2 


1.68 


0.43 


2.35 


2.4 X 


10" 


1.271 


1.8 


X 


1013 


-207 ± 112 


409 ± 121 


1.4 X 


10" 


4.9 X 10" 


-234+ 67 


371+ 82 


-22 




1609-01 


3* 





0.25 


7.82 


8.0 X 


10" 


1.128 


1.30 


< 1013 


-319 ± 114 


407 ± 143 


1.7 X 


1012 


2.1 X 10i2 


-308+ 37 


338+ 58 


6 




1609-01 


4 


-1.68 






















1.4 X 


1012 


1.7 X 1012 


-230+ 88 


409+ 319 


87 




1609-01 


5 


-4.13 






















1.7 X 


1012 


9.1 X 10" 


(-386) 


210+ 94 


10 


00 
CO 


1720-00 


1 


1.63 


0.58 


5.59 


5.7 X 


1013 


1.157 


4.1 


X 


10" 


-451 ± 5 


213 ± 8 


1.7 X 


10" 


6.7 X 10" 


-450+ 4 


212+ 7 


-121 


1720-00 


2* 





0.70 


7.78 


7.9 X 


1013 


1.135 


8.0 


X 


10" 


-397 ± 6 


352 ± 12 


4.4 X 


10" 


1.0 X 1015 


-402+ 5 


342+ 10 


-31 




1720-00 


3 


-2.39 


0.55 


2.90 


3.0 X 


1013 


1.244 


1.8 


X 


10" 


-281 ± 4 


217 ± 7 


1.1 X 


10" 


5.1 X 10" 


-281+ 4 


223+ 7 


110 




1720-00 


4 


-4.03 


0.61 


10.4 


1.1 X 


1014 


1.120 


1.2 


X 


10" 


-280 ± 60 


533 ± 73 


6.1 X 


10" 


1.3 X 10" 


-290+ 43 


565+ 58 


79 




1720-00 


5 


-5.52 


0.77 


9.00 


8.9 X 


10" 


1.128 


1.2 


X 


10" 


-188 ± 12 


407 ± 24 


1.9 X 


10" 


1.4 X 10" 


-196+ 11 


412+ 24 


164 




1836+35 


1 


5.66 






















2.4 X 


10" 


1.6 X 10" 


-259+ 299 


664+ 1516 


-91 




1836+35 


2 


1.68 


0.55 


6.01 


5.9 X 


10" 


1.151 


6.0 


X 


10" 


-400 ± 11 


384 ± 25 


1.5 X 


10" 


1.0 X 10" 


-420+ 37 


375+ 40 


63 




1836+35 


3* 





0.62 


3.71 


3.8 X 


1013 


1.196 


2.2 


X 


10" 


-348 ± 56 


407 ± 29 


1.2 X 


10" 


1.0 X lOi-"' 


-347+ 158 


395+ 175 


-200 




1836+35 


4 


-1.68 


0.46 


2.18 


2.8 X 


1013 


1.240 


7.4 


X 


1012 


-280 ± 51 


202 ± 46 


6.9 X 


1012 


7.0 X 10" 


-259+ 70 


244+ 43 


-183 




1836+35 


5 


-4.59 






















1.8 X 


1012 


1.4 X 10" 


-342+ 294 


347+ 282 


-36 




1836+35 


lb 


5.66 






















1.6 X 


1012 


8.4 X 10" 


34+ 15 


(89) 






1836+35 


2b 


1.68 


0.55 


2.16 


2.4 X 


1013 


1.272 


2.3 


X 


1012 


-135 ± 33 


(89) 


5.8 X 


1012 


3.1 X 10" 


-315+ 65 


345+61 






1836+35 


3b* 





0.62 


2.915 


3.3 X 


10" 


1.226 


2.3 


X 


10" 


-120 ± 7 


226 ± 15 


3.5 X 


10" 


5.6 X 10" 


-121+ 24 


222+6 






1836+35 


4b 


-1.68 


0.46 


4.60 


4.1 X 


1013 


1.185 


2.4 


X 


10" 


-82 ± 21 


242 ± 11 


1.2 X 


10" 


4.0 X 10" 


-71+ 23 


227+18 






1836+35 


5b 


-4.59 






















2.2 X 


1012 


1.2 X 10" 


-10+ 23 


181+ 43 






1929-04 


1 


2.0 


0.33 


2.59 


2.6 X 


10" 


1.259 


3.6 


X 


1012 


-389 ± 24 


191 ± 30 


3.3 X 


1012 


1.8 X 10" 


-391+ 25 


192+ 31 


-179 




1929-04 


2 


1.0 


0.38 


9.19 


9.4 X 


10" 


1.123 


3.1 


X 


10" 


-411 ± 38 


441 ± 55 


1.8 X 


1013 


6.7 X 10" 


-405+ 42 


456+ 65 


-126 




1929-04 


3* 


-0.5 


0.25 


4.39 


4.5 X 


10" 


1.177 


1.8 


X 


10" 


-403 ± 28 


321 ± 50 


8.5 X 


1012 


3.1 X 10" 


-399+ 31 


352+ 68 


-57 




2008-03 


1 


1.87 




2.06 


2.1 X 


10" 


1.314 


1.0 


X 


10" 


-360 ± 12 


193 ± 14 


3.2 X 


10" 


8.5 X 10" 


-225+ 68 


575+ 79 


50 
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Complete text and figures available in a single document at: http://www.physics.ucsb.edu/cmartin/pul 
Figures la-o and Figure 11 available as GIF files on astro-ph. 



Table 3 — Continued 



IRAS Ap Offset Cf T N(Nal) 

(") (A5896) (cm-2) 



(1) (2) (3) (4) 



(5) 



(6) 



DR 



(7) 



N(Nal) V AV min N(Nal) max N(Nal) V AV Vhc 

(cm^^) (km s^-*-) (km s^^) (cm~^) (cm^^) (km s^-"-) (km s^-*-) (km s^-"-) 



(8) 



(9) 



(10) 



(11) 



(12) 



(13) 



(14) 



(15) 



2008-03 


2* 


0.18 


0.52 


9.95 


1.0 X 10" 


1.120 


1.2 X 10" 


-419 ± 


651 ± 21 


2.7 X 10" 


1.1 X 10" 


-398± 44 


680+ 60 


-31 


2008-03 


3 


-1.49 


0.53 


9.50 


9.7 X 10^3 


1.121 


9.8 X 10" 


-391 ± 15 


651 ± 39 


1.1 X 10" 


1.0 X 10" 


-396± 14 


649+ 39 


-126 


2008-03 


4 


-3.18 


0.42 


7.42 


7.6 X lO^-^ 


1.135 


4.8 X 10" 


-404 ± 29 


526 ± 75 


3.5 X 10" 


6.1 X 10" 


-411± 29 


525+ 74 


-82 


2336+36 


1 


3.56 


0.10 


2.45 


2.5 X lO^-^ 


1.277 








1.6 X 1012 


8.6 X 10" 


(-396) 


(328) 


-41 


2336+36 


2 


1.92 


0.22 


3.56 


3.6 X 10i3 


1.210 


9.9 X 10^2 


-478 ± 63 


(328) 


2.3 X 10^2 


3.1 X 10" 


(-396) 


284+ 61 


23 


2336+36 


3* 


0.28 


0.28 


7.44 


7.6 X 10" 


1.136 


2.6 X 10" 


-417 ± 23 


295 ± 28 


2.0 X 10" 


3.2 X 10" 


-413± 26 


299+ 30 


-19 


2336+36 


4 


-1.36 


0.43 


3.04 


3.1 X 10" 


1.234 


1.2 X 10" 


-300 ± 158 


460 ± 269 


9.6 X 10" 


2.6 X 10" 


-381± 56 


303+ 66 


-115 


2336+36 


5 


-2.99 


0.22 


2.15 


2.2 X 10" 


1.304 


6.2 X 10^2 


-370 ± 191 


458 ± 344 


3.6 X 10" 


1.8 X 10" 


-345+ 239 


376+ 438 


-19 


2336+36 


6 


-6.57 
















2.5 X 10" 


1.3 X 10" 


-747+ 1081 


399+ 3408 





Note. — (1) Galaxy; (2) Aperture number from Figure 1, where asterisk denotes aperture nearest the R-band nucleus; (3) Offset of aperture from nucleus; (4) Estimated 
covering factor; (5) Fitted optical depth of the dynamic component; (6) Column density of the dynamic component from optical depth fit in col (5); (7) Doublet ratio 
chosen for intensity fit based on value in column (6); (8) Column density of the dynamic component from intensity fit with doublet ratio Dn from column (7); (9) Doppler 
shift of the dynamic component (relative to CO velocity) from the intensity fit with Dji from column (7); (lO)FWHM of the dynamic component from the intensity fit 
with Dfi from column (7); (11) Maximum column density from the intensity fit with variable D/j, where the maximum of columns (11), (6), and (8) should be taken as 
the best value; (12) Minimum column density from the intensity fit with variable Dn, the minimum of columns (12), (6), and (8) should be takend as the best value; 
(13) V from the intensity fit with variable D^; (14) FWHM from the intensity fit with variable D^. Throughout this table, parameters listed in parentheses were held 
constant to get an adequate fit. 



hJ^ 
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Fig. 1. — {Top Left) Position of the apertures on an R band image from Murphy etaL 
(1996). North is up, and east is to the left. The sht is 20" long. {Top Right) The Ho; plus 
[Nil] AA6548, 6584 emission lines along the slit; wavelength increases to the right. Spectral 
apertures and slit orientation are indicated. {Bottom Right ) Extracted Ha spectra with 
apertures and slit orientation marked. {Bottom Left) Extracted spectra at Na IAA5890, 96. 
The nuclear spectrum (red), aperture number, and slit orientation are marked. The blue 
line indicates the systemic velocity of the A5896 line. Two Na I doublets are plotted ; the 
velocity of one is fixed at the systemic velocity (dotted), and the velocity of the dynamic 
component is fitted (solid). The line shapes of two galaxies required a minimum of two 
dynamic components, and the second dynamic component is shown by a dashed doublet. 
The emission line is He I A5876. (a) IRAS10565+2448 - Blueshifted, interstellar absorption 
is detected everywhere continuum is detected. In contrast, the Ha emission is redshifted 
(blueshifted) west (east) of the nucleus consistent with galactic rotation. The companion 
to the southeast shows no Ha emission (i.e. little star formation) but is apparently covered 
by the outflow. The Na I linewidth and covering factor decrease away from the nucleus, 
but the average outflow velocity remains near 300 km s~^ . (b) IRAS 17208-0014 - The Na I 
absorption is highly blueshifted across both of the merging galaxies; the outflow speed is 
400 km s~^ at the nucleus. The Ha emission reveals rotation across the main galaxy and 
the prograde nature of the encounter. In contrast to the IRAS10565+2448 spectrum, the Na I 
absorption shows the same velocity gradient as the Ha emission here. Additional galaxies 
are shown in the electronic edition as follows: (c) IRAS00153+5454, (d) IRAS00188-0856, 
(e) IRAS00262+4251, (f) IRAS03158+4227, (g) IRAS08030+5243, (h) IRAS10494+4424, 
(i) IRAS11506+1331, (j) IRAS15245+1019, (k) IRAS16090-0139, (1) IRAS16487+5447, (m) 
IRAS18368+3549, (n) IRAS19297-0406, (o) IRAS20087-0308, (p) IRAS23365+3604. 
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FIGURE 2a-o available as f2a-o.gif 

Fig. 2. — Gas velocities measured across the ULIGs in Ha (circles) and Na I. (squares). 
Apertures are marked in the top left panel of Figure 1. The Ha emission traces the pro- 
jected rotation of the merging galaxy pair along the slit. The Na I absorption traces the 
speed of the cool outflow. Aperture numbers, slit orientations, and the location of double 
galaxies are marked for comparison to Figure 1. The galaxies are: (a) IRAS10565+2448, 
(b) IRAS17208-0014, (c) IRAS00153+5454, (d) IRAS00188-0856, and (e) IRAS00262+4251, 
(f) IRAS03158+4227, (g) IRAS08030+5243, (h) IRAS 10494+4424, (i) IRAS11506+1331, 
(j) IRAS15245+1019, (k) IRAS16090-0139, (1) IRAS18368+3549, (m) IRAS19297-0406, (n) 
IRAS20087-0308, (o) IRAS23365+3604. Note that two dynamic absorption components were 
fitted in IRAS10565+2448 and IRAS18368+3549 as discussed in Martin (2005), so there are 
two Na I measurements at each position. Uncertainties in the Ha points are described in 
S3. 2.1. 
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AVAILABLE AS f3.gif 

Fig. 3. — Comparison of Na I and Ha velocity gradients along the slit. The total Na I 
velocity gradient is compared to the Ha velocity gradient: (Left) across the entire system; 
(Center) across one galaxy - i.e. galactic rotation; (Right) between the nuclei, i.e. orbital 
motion. The small solid squares denote single nuclei systems. The velocity shear across the 
outflowing gas is comparable to that of the bulk interstellar medium, which reflects galactic 
and orbital rotation. 
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AVAILABLE AS f4.gif 

Fig. 4. — Position-velocity diagrams for Na I absorption grouped by velocity gradient, (a) 
The absorption velocity is nearly constant along the slit. If the outflow diverges from the 
central starburst region, these gradients are consistent with projection along the major axis of 
the gaseous disk, (b) The absorption velocity jumps near the starburst nucleus. The conical 
wind model shows this type of shear when projected along the minor axis of the gaseous 
disk, (c) The velocity of the absorbing gas presents rotation. Rotation is not consistent with 
a flow diverging from a nuclear starburst. A larger launch region and/or additional torques 
are indicated. 
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FIGURE 5 available as f5.gif 

Fig. 5. — Masses of the cool outflows estimated by integrating the Na I column density 
over each galaxy. The non-linear scaling of the Na Icolumn density with the line equivalent 
width creates a large uncertainty in any one mass estimate. The columns were estimated 
from the apparent optical depth and could be much higher if strongly saturated components 
are present. The conversion from N(Na I) to total gas column is assumed to be a factor 
of ~ 10 (see text for details). The results suggest > 10® Moof cool gas is carried by ULIG 
winds. The small range of far-infrared luminosity indicates the baseline in star formation 
rates among the systems shown is small. The symbol denotes the dynamical age of the 
merger: double nuclei (solid), single nuclei with prominent tidal tails (open), and single 
nuclei (crosses). 
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FIGURE 6 available as f6.gif 

Fig. 6. — Mass of the cool outflow vs starburst dynamical mass. Outflow mass is from 
Table 1. Dynamical mass estimated from CO linewidths by Solomon etal. (1997). Triangles 
denote wind masses from Schwartz & Martin (2004) and Schwartz (2005), molecular gas 
masses for M 82 (Weiss etal. 2001) and NGC 1614 (Scoville etal. 1991), atomic gas mass 
from Kobulnicky etal. (1999) for He2-10, and the atomic gas masses tabulated in Schwartz 
& Martin for the remaining systems. Taking the dynamical mass in the starburst region 
as a measure of the total gas mass, the diagonal lines show outflow gas fractions of 100% 
(dotted), 10% (dashed), and 1% (solid). 
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FIGURE 7 available as f7.gif 

Fig. 7. — Color excess versus Na I equivalent width. Uncertainties are dominated by correc- 
tions for stellar absorption; and the full range allowed is shown by the error bars. ULIGs 
with poorly constrained reddening are shown by open squares. Filled circles are LIGs from 
Heckman et al. (2000). The systems with stronger interstellar Na I tend to be more reddened. 
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FIGURE 8 available as f8.gif 

Fig. 8. — Covering fraction of absorbers vs. logarithmic extinction at H/9 . Data are nuclear 
measurements for the ULIGs (squares) and the LIRGs from Heckman etal. (2000) (circles). 
Galaxies with poorly constrained reddening are shown by open symbols. Systems with nearly 
black absorption lines are highly reddened. 
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AVAILABLE AS f9.gif 

Fig. 9. — Estimated kinetic energy of the cool outflow vs feedback energy supplied by 
10 Myr starburst. The lines denote 100% (dotted), 10% (dashed), and 1% (solid) of the 
energy supplied by a 10 Myr old starburst. Triangles denote dwarf starbursts from Schwartz 
& Martin (2004), who used a larger corretion for sodium ionization and depletion than 
assumed here for the ULIGs (see text for details). 
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FIGURE 10 available as flO.gif 

Fig. 10. — Mass loss efficiency in cool phase vs terminal outflow speed. Estimates of the mass 
flux are normalized by the star formation rate (of 1 to 100 M0 stars), and the velocities of 
the most blue-shifted absorption are estimated from Paper I for ULIGs (square points). The 
supershell mass from Schwartz & Martin (2004) was divided by a dynamical age of 10 Myr 
to estimate mass flux (triangles). The dotted lines illustrate models for energy-driven winds 
in which 100, 10, or 1% of the mechanical power is thermalized. Dashed lines describe the 
shallower relations for momentum-driven winds in which 100, 10, or 1% of the mechanical 
momentum goes into the wind. The available mechanical power and momentum (per unit 
SFR) is taken from a Starburst 99 population synthesis model with continuous star formation 
for 10 Myr, solar metallicity, a Salpeter initial mass function, and lower (upper) mass limits 
of 1 Mq (100 Mq). The ULIG data points (squares) are consistent with outflows driven 
by a few percent of the mechanical energy or ~ 10% of the momentum. The variation in 
mass-loss efficiency with terminal velocity distinguish these models; but this interpretation 
is subject to how the dwarf starburst points are estimated. 
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AVAILABLE AS fll.gif 

Fig. 11. — Geometrical toy model of a conical wind leaving a nuclear starburst region of 
radius Rq. The rotational velocity of the galactic disk, Vc, and the radial outflow velocity, 
Vw, are projected onto the observer's detector at upper right. The observer is positioned at 
infinity in the yz-plane. The disk is inclined by an angle t with respect to the observer. The 
red sightline intersects the major axis of the disk, which is aligned with the x-axis in the 
diagram. 
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FIGURE 12a,b available as fl2a,b.gif 

Fig. 12. — Projected (rotation plus outflow) velocity of toy model for constant-L, wind (red) 
and constant-l), wind (blue). The radial wind speed is taken as twice the rotation speed of 
the disk; both are constants. The half-opening angle of the wind cone is 60° . The thick-solid, 
solid, and dashed lines denote viewing angles of 55° , 35° , and 15° , respectively. There 
is no wind inside the launch radius, Rq. (a) Galaxy major axis: The step function shows 
the rotation curve of a disk rotating at Vc- Points mark the starburst ring at Ro, where the 
wind is launched (green). Compared to the disk, the wind presents little velocity gradient 
regardless of viewing angle, (b) Galaxy minor axis: No rotation is observed along the minor 
axis. The wind cone produces a velocity jump because sightlines on the near-side of the disk 
are more aligned with the outflow than those on the far side, which are nearly perpendicular 
to it (red lines). 
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